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equations, is numerically complex and time consuming. Also 
as tlie Reynolds num][>er Increases, It l>ecomes increasingly 
difficult to capture the thin shock through the merged 
layer computation. One of the main aim of the present 
Investigation Is to increase the domain of the validity of 
the solution and to avoid the numerical Instabilities 
associated with the Integration through the shock. To 
achieve this, a shock layer model with full Navler-Stokes 
equations is proposed In this work. The domain of validity 
Is Increased because It retains all the essential terms in 
the equations. The numerical Instabilities do not limit the 
solution because the Integration through the shock 
transition zone Is avoided by employing a shock layer 
model. 

With the emergence of real gas effects at high 
temperatures, the problem of simulating the flow-field 
becomes more complex by any of these approaches due to the 
presence of mass production terms in the species 
equations. Hon- linear dependence of the production term on 
temperature and concentrations introduces stiffness Into 
the equations, which usually leads to non-convergence and 
failure of the numerical computations. Sometimes It also 
leads to a violation of the overall mass conservation. 
These problems become more severe as the Reynolds number 

increases. Hie present investigation also ainm at studying 

*• 

these limitations and removing them to the extent possible 
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l>y a new shocK layer formulation wim full Navler-StoKes 
equations. 

In the present investigation, the full Navier-StoKes 
equations are used to describe the hypersonic flow near 
the stagnation point of a blunt body both for perfect and 
chemically reacting gases. The equations are solved 
between the body and the shock. The surface slip 
conditions for the multi> component gas mixture are used at 
the wall. The mass production rates within, the shock 
transition zone are neglected but species diffusion at the 
shock is accounted for in the shock conditions. Modified 
Rankine-Hugoniot relations are used as boundary conditions 
at the shock. Using the concept of local similarity, the 
governing equations are specialized for stagnation line 
flow. The resulting system of non-linear coupled ordinary 
differential equations with split boundary conditions is 
solved numerically using a special finite difference 
scheme called Unified Difference Representation. The 
equations are solved both for perfect and chemically 
reacting gases. 

It is demonstrated that the flow field simulations 
using the full Navler-Stokes equations in the stagnation 
zone of a blunt body provide uniformly accurate results 
in low as well as high Reynolds number regime. At low 
Reynolds numbers it shows a, good agreement with the 
experimental data. The computed results reproduce 
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experimentally observed featiires of low density flows 
llKe, Increase in impact pressure and recovery temperature 
at stagnation point wltn decrease in Reynolds number, 
wblcli tbe earlier approximate methods fall to predict even 
qualitatively. The computed beat transfer coefficient 
shows good agreement with boundary layer results at higher 
Reynolds numbers. The predicted profiles agree well with 
the available experimental measurements for density, 
pressure and temperature profiles across the shock layer 
for the low density regime. 

The analysis shows that the wall temperature. In 
addition to the Cheng's rarefaction parameter, plays a 
vital role in determining the heat flux to the wall.’ For 
chemically reacting flows In addition to the Reynolds 
number, the free stream velocity also influences the heat 
transfer coefficient to the wall and needs to be accounted 
for explicitly. The concentration profiles are in good 
agreement with the merged layer predictions. This good 
agreement holds even at low Reynolds numbers where the 
shock wave is thick. The special treaitment of mass 

production terms along with full Kavler-Stokes equations 

\ 

In the shock layer has resulted in the Increased range of 
Reynolds numbers for which non> equilibrium computations 
converge. These factors establish the superiority of 
present shock layer formulation with the full Havler- 
Stokes equations over the current viscous shock layer 







1 NTROIDUCO? I ON 


1 IMHOPOCTIOl 

nie problMB of bypersonlc flow past a blunt 
azlSTometrlc body bas attracted considerable attention due 
to Its application to the design of re-entry vehicles. 
Usually, aerodynamic drag Is used to slow down the 
vehicles entering the atmosphere with hypersonic speeds. 
During the atmospheric entry phase of a vehicle at 
hypersonic velocity, considerable work Is done by the 
vehicle on the air surrounding it due to viscous effects. 
This work done by the vehicle on the surrounding 
atmosphere exhibits itself as thermal energy and results 
in increased temperature of the air surrounding the body. 
A part of this heat Is transferred to the vehicle Itself 
and increases the temperature of the body. In order that 
the vehicle survives this thermal environment, the body 
has to be protected against this heat load, which Is quite 
Intense. This protection Is required from two 
considerations — (l) too much heating is likely to change 
the aerodynamic shape due to ab 1 at 1 on/me 1 ting of the 
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veli.lcle's outer skin resulting in chanpeci aerodynamtc 
cliaracteri sties and performance of tfce ver^icie fcre « 
tliis increased temperature of tHe vehicle skin »fiil 
in change in the environment of the electronic equipj.jcri-'.s 
and effect their proper functioning. Thus it is essential 
that one has a good estimate of heat flux load coming on 
the vehicle. so as to he ahle to design an efficient 
protection system against this aerodynamic heating ar.dl 
ensure desired vehicle performance. 

As vehicle descends in the atmosphere, the 
surrounding air changes its characteristics from the free- 
molecular flow regime in the upper atmosphere to the 
continuum regime in the lower atmosphere. The free- 
molecular regime is well investigated from kinetic theory 
considerations. The other extreme of the highly dense air, 
where shock wave is extremely thin and the viscous effects 
are confined to a narrow region near the hody surface, 
lends itself to the classical boundary layer treatment. 
Between these two limits, the situation exists in which 
air is sufficiently rarefied resulting in hreakdowii cf 
usual boundary layer approximations but at the same time 
mean-free path is not large enough for the free-molecuiar 
concepts to apply. This flow regime is referred to as the 
intermediate regime. The main features of this regime ai-e 
the thickening of the boundary layer and the shock wave 
and their subsequent merger. Usually this regim: if 



furtiier stilxli victed i>aseci on tiia differer.t tr^ecr 
models conslcered necessary xo CLScribe tnc flow. Vt-rco-cs 
authors have subdivided tins intern^c diate re* irr-o- 

differently depending ujon tlie flow phenorTit-nc-n anJ. tlie 
range of validity of the reEime. Kays and Probstein (iOCt] 
divided the intermediate regime into the foliov int 

regimes : 

(I) First order collision regime 

(II) Transitional regime 

(ill) Fully merged layer regime 

(IV) Incipient merged layer regime 

(V) Viscous layer regime 

(VI) Vorticity interaction regime 

However, Ahouse and Bogdonof f ( 1969 ) divided this 
intermediate regime into following two regimes: 

(I) Transitional regime 

(II) Merged layer regime 

The transitional regime defined hy Ahouse and 
Bogdonof f includes first order collision regime and 
transitional regime as defined hy Hays and PiohS'-ein, 
while merged layei' regime of Ahouse and Bcgcoaoff 
encompasses the fully merged layer regime, incipient 
merged layer regime, viscous layer' regime aad vorxicity 
interaction regime of Hays and Prohstein* In this verh the 
classification given hy Ahouse ana Bogaonoff is fol loved. 
The thicKening of the boundary layer ever the body 



surface is associated with decreasing Reynolds number «nd 
increasing Kach number. Increased Hach number results in 
increased temperature in the flow-field. This increased 
temperature leads to higher viscosity resulting in 
flow-field being dominated by viscous effects. At the same 
time, increased Hach number results in decreased shock 
stand-off distance. Decreasing Reynolds number is 
associated with increase in viscous forces in comparison 
to the inertia forces. Ihis increased viscous effect lead 
to increased boundary layer thickness. Thus at higher Hach 
number and lower Reynolds number regime, the entire 
flow-field is dominated by viscous effects and there is no 
inviscid flow region between the body and shock wave. 

Another important feature of hypersonic flow is the 
emergence of the real gas effects, which means that the 
thermodynamic properties of the real air deviate from 
those of an ideal gas with constant specific heats. These 
deviations result from various physical processes 
occurring in uhe air at high temperature, namely, 
eircitation of molecular vibration, dissociation ionicction 
ana chemical reactions. These processes occur at a finite 
raie, iRhich can be characterized by relaxation time or 
relaxation lengths in a flow of a given velocity. If all 
the relevant relaxation lengths are very much smaller than 
the smallest flow field dimension of interest, flow may be 
regarded as being xn thermo -chemical equi 1 ibri\im. If 



relaxation lengtl^is and flow dimensions are comparable, 
departure from tlxermo- chemical equilibrium will occur. The 
resulting non- equilibrium flow of real gas may influence 
the aerodynamic cheu'acteri sties. The analysis of non- 
equilibrium air flow requires detailed Information on the 
mechanism and the rate of high temperatiire processes. If 
all the relevant relaxation lengths are much greater than 
the largest flow dimension, the flow may be regarded as 
frozen flow and chemical composition remain constant. 
These relaxation lengths are very large at higher 
altitudes and quite small at lower altitudes. Thus during 
re-entry of a vehicle, non- equilibrium flow must occur. 

From continuum considerations, the existing 
theoretical analyses can be divided into three main 
categories. Firstly, there is the higher order boundary 
layer theory, developed systematically from the theory of 
singular perturbations. The classical boundary layer 
theory of Prandtl is the first approximation of this 
theory. Higher order terms provide corrections to the 
first approximation. This higher order theory is valid 
only If there is small departure from the classical 
boundary layer limit. The second approach to the blunt 
body problem is the two thin layer theory, where the flow- 
field is divided into two zones-(i) shock transition zone 
< shdcK structure) and (li) shock layer that is the region 
between the shock and the body. In this approach, the 
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mass, momentum and energy fluxes tangential to tiie shock 
are assumed to he negligible in the shock transition zone. 
Equations similar to Prandtl's boundary layer equations 

are used to describe the flow field in the shock layer. 

( 

These equations are uniformly valid for both inviscid flow 
and viscous flow. Transport effects just behind the shock 
wave - are taken into account and consequently the classical 
Rankine-Hugoniot shock conditions are modified to account 
for the transport effects. These conditions are derived by 
integrating the ordinary differential equations valid in 
the shock transition zone. Investigations of the flow 
field by this model are -found to be successful in 
predicting surface quantities such as heat transfer rate 
well-* into the transition region, but fail to give correct 
behaviour of surface pressure and recovery^ factor at low 
Rpynnid's number regime. It also fal!~ *•' d<»«cr.t>''“ the 

detailed flow-field accurately at low Reynolds numbers. It 
is surmised in this work that this failure to predict the 
correct behaviour of surface quantities at lower Reynolds 
number regime is due to the simplification carried out in 
the: Navier-Stokes equations and not due to failure of 
shock layer model as such. Prom continuum theories, the 
last: approach to the problem is the direct numerical 
integration of Havier-Stokes equations from the surface to 
the free-stream. This method is numerically complex and 
time consuming. Also as the Reynolds number increases, it 
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becorsfes i^^’^t'easingly difficvilt to caj ture tiie tiiiri &hocii 
hj direct r Integration. Several attempts have been made to 
investigate non-equllibri-um flow using tbese models. Tbe 

ft 

governing e'^uations of tbe non*equi librium flow are stiff 
due to mass production terms in tbe species continuity 
equations. Tbis stiffness leads to convergence problem and 
sometimes to violation of overall mass conservation. 
Currently Honte Carlo simulation is being widely used to 
simulate tbe flow field in tbe transition regime. 

Tbe main aim of tbe present dissertation Is to 
improve the continuum solutions in tbe stagnation region of 
a blunt body under bypersonlc conditions witb and without 
tbe real gas effects. Improvements are required to 
increase tbe range of validity of solutions and to remove 
tbe instabilities associated witb tbe tbin sbocR. For 
cbemically reacting flows improved method is required to 
deal witb stiffness of tbe species continuity equations 
and to improve tbe convergence. Real gas effects include 
excitation of various degrees of freedom ( exhibited as 
variation of specific beat with temperature ) and chemical 
reactions. Surface slip and shocK slip conditions are 
considered to account for rarefaction effects. 

2 Aaaessiaept of Existing Theories and Literature 

Existing continuum approaches to tbe hypersonic blunt 

; 

body problem can be divided into three main categories. 
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4. ' f- ( I c 

( 1 ) HifJier oraei-' bo-a.,icif.ry iayt.-r xheory 

Ui) Two tliin layer theory/ £;iOCit layer theoi-y 

(iii) Merged layer theory 

However, currently Kinetic theory is being used to 
prf'dict the flow-field well into transition regime. This 
approach is termed as 

(iv) Direct Simulation Monte Carlo (DSHC) technique 

A survey of continuum methods was made by Cheng 
(1966) and Jain (1966). A critical evaluation of the work 
done along these approaches, leading to the motivation of 
the present investigation, is given below. 

Higher Order Boundary Laver Theory 

The classical boxindary layer theory of Prandtl is 
embedded as the leading approximation In a systematic 
asjnmptotic expansion of the Havier-StoKes equations in an 
inverse power of the Reynolds number. The method of inner 
and outer expansions, also called the matched asymptotic 
expansion ( Van Dyke 1964a ) is used here. Different types 
of expansions are used in different regions of the flow, 
depending on the order of magnitude of the thickness of 
the regions. So far, calculations are made up to third 
order boundary Ifiyer, wl^h Prandtl 's boundary layer as tlve 
leading first approximation. ExT-austive work on the second 
orcer boundary layer is done by Van Dyke (1962, 1963, 

1964b), Rott and Lenard (1962) and Lenard (1962) among 
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others. CersF '..f.tioris Vrrf c&rritia ott for third order 
hounc.ery l£>e.:- hy ilao (19&*»a) for cold va; 1 a:id by Oberai 
11964'; for an insulated wall. 

The procedure followed in the higher order boundary 
layer theory is as follows. The flow quantities are 
expanded asymptotically for large Reynolds number in the 
shocK wave, the Inviscid region and the boundary layer. 
The variables are stretched appropriately for the 
different flow regions depending upon the thicKness of the 
region. For proper stretching of the coordinates and 
variables, estimates of the thicKness for the various flow 
regions already avail2d>le from previous studies, are used. 
These expansions are substituted in the Havier-StoKes 
equations and the coefficients of like powers of the 
stretching parameter are collected resulting in different 
sets of equations valid in different regions of the flow. 
These equations are solved according to the following 
scheme: 

(i) First Order Boundary Layer Approximation: - (a) first 
the basic inviscid flow is computed and (b) the first 
order boundary layer equations are solved matching with 
the basic inviscid flow. 

(ii) Second Order Boundary Layer Approximation:- ia) The 
second order equations in the inviscid region are solved 
takihg into account the displacement thickness caused by 
tile , order! boundary layer, (b) the second order 



; c 

t-: > iaisr equaticnr fc.:-e tolyt 1 ir, itching w^th thfe 
i !.Vi SClh flow. 

iaix) Third, Order Boundary Layer Approximation*. - (a) The 
fxr&t order equations for the shocH structure are solved 
matching with the free-stream ahead and basic invlscid 
flow behind, (b) the third order equations in the inviscid 
region are solved talking into account the modified shock 
conditions, which depend on the first order shock- 
structure solution and the shock curvature, (c) the third 
order boundary layer e<xuations are solved matching with 
the third order inviscid flow, and finally, (d) the third 
order shock structure equations are solved. 

% 

For any given Reynolds number, the afore mentioned 
solutions for the various components can be combined to 
give a composite solution that is uniformly valid 
throughout the flow field. 

The higher order boundary layer theory can describe 
the flow accurately only if the higher order terms are 
small in comparison to the basic Inviscid and first order 
solutions. This implies that Reynolds number is moderately 
large and interaction between the boundary layer, inviscid 
flow and shock structure is not significant. This analysis 
allows the identification of several second order effects 
rjid calculation of their influence on the flow field 
individually. Van Dyke U96£) identified seven second 
order effects, namely, longitudinal curvature, transverse 
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CiityK',':, €. 'vfiri.al eraa^ent of et tJiSipy, slip velocity ar.d 
f J'etnr-e at tlie w&li. TiAe exteriial vorticity is 

relaxed to th.e external entropy gradient and enthalpy 
gradient tlm otgh the Crocco's relation. All these effects, 
except for velocity slip and tenqperature jump, are not 
independent of each other in strict sense and thus it is 
better to evaluate them concurrently. Van DyKe (1963) 
evaluated these Influences for stagnation point heat 
transfer rate for a sphere at H® = 4 for perfect gas with 
V = 7/5, Prandtl number 0. 7, viscosity proportional to 

temperature and surface temperatxire being O. 5 times 
stagnation temperature. He gave the following expression 


for 

heat transfer 

rate: 
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Similar expression is available for the case when 
surfece temperature is 0.2 tiroes stagnation temperattire. 
Here EOg is the Reynolds number based on the post-shocK 
conditions and nose radius. It is seen here that the 
external vorticity has major influence on the stagnation 
joint heav flux, a detailed discussion on the various 
second ord. r effects can be found in the review papers by 
Cheng (1966), Jain (1966) and Van Dyke (1969), 



T... M-cc: bouridai'i' icy^r tn-rcrv dc r.ot ^::ve 

.'vT-uit r-6-vits fcr lew E;«.y;'iC'ACii i’lew, ^ 

fc;r-ire flow fisld isetween the body and the sbocK is 
inf Itt'jnced by vise '•us effects. The boundary layer tbeory 
en’v'isugfs an essentially inviscid zone between tbe sbocK 
wave and tbe boundary layer, even thougb in tbe tbird 
order boundary layer theory viscous terns appear in this 
zone. In this flow regime, the departure of flow 
quaitities from the classical boundary , 1 ayer .solutions is 
large and thus cannot be treated as a perturbation. 
Koreover, the method of treating the various zones in a 
piecewise manner is no more valid in this flow regime. 
Indeed, Kao (1964b) noted that for Reg less than about 
100, the heat transfer rates predicted by the third order 
boundary layer theory deviate appreciably from the 
experimental results. The failure of the higher order 
boundary layer to describe the low Reynolds nuna>er flow 
exi-resses itself in interesting ways. Anomalous situations 
such as the boundary layer overtaking the position of the 
sLock or the thick shock wave encompassing even the body 
(Kso IS 64b) occur during the computation of low Reynolds 
nuyher flow. In such cases, it is difficult to construct 
£.:.y meaningful com:, ©site solution. 

For the case of real gas and chemically reacting 
slow, most of the studies are confined to the classical 
bev-ndery layer analysis and hi r her order boundary layer 



't-uS r.r>t been ct-'Ve loj-e i. 

Fay &it' R.'ddfril (l?5t) analys&e the influence f»* 
cbcr-ically reacting flow on tlje stagnation point beat fius 
ana came sut vitb a correlation for surface be ,t transfer 
rate for frosen and equilibrium flow. This analysis is 
valid for large Reynolds ntimber, when the viscous effects 
are confined to the region neeir the wall and an inviscid 
flow region separates the boundary layer and the shocK. 
The analysis is based on the assumption that air is a 
binary gas and that the Prandtl and Lewis numbers are 
constant across the boundary layer, though the 
concentration changes across the boundary layer due to 
temperature variation . This correlation expresses the 
stagnation point heat flux in two parts — (l}heat flux due 
to perfect gas and (ii) contribution due to species 
diffusion both for frozen and equilibrium flows. Blottner 
(1964) analysed the non-equilibrium boundary layer for a 
binary gas consisting of atoms and mol^ cules. This 
analysis required the inviscid non- equilibrium flow 

conditions on the surface as input. Ihe governing 
equations were solved numerically using a marching 
procedure starting from stagnation point/leading edge of 
vhe body. Tlie analysis shows that in the non- equilibrium 
boundary layer the maximum temperature could be nearly 4. 5 
tires higher than the corresponding temperature for the 
equilibrium boundary layer. 
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Tr.c? Biixn f eatwi e of t:iis l]rieory is to <Jiv,<5e tJie 
eritire flov field into two asijoining regions, namely the 
shocK layer and the shock wave < or the shock transition 
zone), hoth of which are assucied to he thin. In the shock 
transition zone, fluxes of mass, moQientum and energy 
parallel to the shock are neglected. This assumption, 
together with the thin layer approximation, reduces the 
governing equations in the shock transition zone to a set 
of ordinary differential equations similar to that for a 
plane shock. The governing equations for the shock layer, 
that Is, the region between the Inner edge of the shock 
wave and the body, embody differential equations of botl 
the inviscld flow and the classical boundary layer theory. 
These are obtained from Havier-Stokes equations by using 
the thin layer approrimatlon and a high shock compresslor 
j-atio assumption by retaining the terms of the order oi 
e = (y-i)/£y. These equations are parabolic in nature 

along the length of the body and their mathematical 
analysis is similar to the problem of classical boundars 
layer theory. 

The features of t,he viscous shock layer theory are 
similar to two thin layer theory. However, in this case 
€quc-.t*cns ai'e solved in the shock layer only, that is the 
Sifdficn between the body and the shock. The equations ir 
11*7 Shock transition region are not solved. While deriving 


T' ■ f f.'i - rn-^r.r ec/jaixcr^£., tJiin ciiocK ls,/er fef rursptiori :s 
r.< ct' cr.c X‘Vi*:r-s < i the or'C.e"" of {i/ke) aj'e rete.irieci in xite 
K£.vie-:'-stoXes equations. These equations are uniforaji}’' 
vaiid in inviscid. flow region and viscous flow region near 
the wail. 

An essential feature of the two thin layer model is 
to preserve the viscous and energy transport effects Just 
behind the shock. Ihe higher order boundary layer theory 
failed to do it. The shock transition zone equations, when 
integrated across the shock, lead to modified form of 
Rankine- Hugoniot relations for the post-shock flow 
variables. According to these modified relations, the 
tangential velocity and the total enthalpy are not 
conserved across the shock but are controlled by the 
transport effects behind the shock, which are not small 
for low Reynolds number flow. These transport effects are 
proportional to the gradient of velocity and temperature. 
Being similar to surface slip conditions, these conditions 

i 

are referred to as shock slip conditions (Probstein and 
Pan i9€£). 

It is interesting to note that, in the two thin layer 
t?'»ory, the flow in the shock layer together with the 
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incsepencSent of the flow in the shock transition zone, 



OrCu if j, t-w KE.rJ'".iri£-Kui oniot cor.citioni eris t'zie appi-opriate 
coi.iitic.js at the ir-ody surface are sujfxcient to solve the 
shod; layer equations completely. However, the flow in the 
;:hoci: iransition zone depends cn the shock layer solution. 
Thi- shock layer solution provides the explicit boundary 
conditions for invegrating the first order differential 
equations governing the shock transition zone. 

With this basic formulation of the two thin layer 
theory, Cheng <1961) investigated analytically the shock 
layer flow in the stagnation region in the merged layer 
regime. It is noticed that the use of shock slip 
conditions alone ( without surface slip conditions ) lead 
to reduction of surface heat flux. Thus surface heat 

I 

transfer rate and skin friction values tend towards the 
free-molecular values of heat transfer and skin friction 
for the case of iinit accommodation coefficient. Further 
Cheng and Chang (1964) evaluated an analytical solution 
for the shock transition zone in the Newtonian (y -> 1 ) 
and the hypersonic ( -> oo ) limit. This solution 
compc-res well with the results of Levinsky and Yoshihara 
(1962.), where direct Integration through the shock was 
caj-ried out for monatomic gas. Davis (1970) solved the 
VISCOUS shock layer equations on different hyperboloids. 
Slip conditions are incorporated as higher order terms in 
a perturbation scheme. The shock slip conditions have been 
used to account for transport effects at the shock. 



tl-ie sliocii s-lip ccrsdxtior.r used by Davis are 
iKCcni:- : 1 witn tixe roverning ecjsiions. Tiiese shoe;: 

siiP »^c.nciitions do no^ account for change in post-shocK 
preicure <iue xo cliange in post-sJiocK ncrml velocity. 

Basil (1964) formulated tlie iiypersonic viscous blunt 
body problem using asimptotic expansion of the variables 

in ( i/Re ), ( i/Ka^) and . Busb divides tbe 

y + 1 

shocK transition zone into three subregions, namely (a) 
the outer region, (b) the middle region, and (c) the inner 
region. Including the shock layer region. Bush analyses 
the flow on the basis of four distinct regions as against 
only two used by Cheng in the two thin layer model. Cheng 
(1966) established that the equations governing the shock 
transition zone in the two layer model are equivalent to 
the comrosite equations of Bush and are uniformly valid in 
the enure shock transition zone. 

Regarding the real gas flow with chemical reactions 
around blunt bodies, most of the work has been carried out 
in the viscous shock layer model instead of the two thin 
Uyer model. Blottner (1969) analysed the viscous shock 
layer flow at the stagnation region of a blunt body for 
ca] iy reacting flow. In this work, air was considered 
consiiting of seven species with seven chemical 
reacuions, Tne analysis considered fully catalytic wall, 
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but did jiot account for eiilitr wail slip or transpo'»'t 
effects at tbe sbocK. For evaluating tbe post-sbocK 
conditions, chemical composition was considered to be 
frozen as in the free stream. T.ie governing equations were 
solved using finite difference scheme. Computations showed 
that this analysis predicts higher ten«>erature at the 
shock than the corresponding ecjui librium flow temperature. 
This temperature decreases as one moves towards the wall 
due to dissociation of gas, but the gas has not reached 
the equilibrium before boundary layer is reached. The 


computations 

also 

show that the assumption 

of 

inviscld 

equilibrium 

flow 

for the boundary layer is 

not 

correct 


even for Reynolds number corresponding to an altitude of 
45.72 km. Kang and Dunn (1972) used shock layer analysis 
for hemisphere -cone using eleven species model with twenty 
six chemical reactions for simulating the real gas flow. 
However, the analysis followed the two thin layer model of 
Cheng and assumed constant Prandtl and Lewis numbers and 
binary diffusion. Hodified Rankine-Hugonlot relations were 
used to account for transport effects behind the shock. 
These served as the initial conditions dPor the shock 
transition zone. The method of integral relations was used 
to obtain the solution of the governing equations. The 
analysis shows that, at down stream stations, the 
temperature profile overshoots the post-shock temperature. 

Moss (1971) extended the analysis of Blottner and 
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analysed the chemically reacting and equilibrium flow 
field around hyperboloid. This viscous shock layer 
analysis considers only 5 species air model and does not 
consider Ionised species and ionisation reaction. The 
analysis includes the effect of multi-component diffusion 
and mass injection. 

Computed results of Moss show that the non- 
equl librium model predicts heat transfer rates which are 
lower than both the frozen and equilibrium predictions. 
Non -equilibrium heat flux predictions for the stagnation 
point are approximately 60X of the corresponding 
prediction for frozen and equilibrium flow for the cases 
considered In the' analysis. Results also show that, for 
heat flux predictions, the use of veuc'lable binary 
diffusion coefficient { atomic oxygen In molecular 
nitrogen } compares favourably with multi -component 
diffusion, though there is some difference In the 
concentration profiles near the wall In the two cases of 
binary diffusion and multi-component diffusion. The 
computations fall to converge for the cases when the post- 
shock Reynolds niimber is 500 and above. 

There have been atteim>ts to extend this analysis of 
chemically reacting flow to three-dimensional bodies. 
Hl.ner and Lewis (1975), Elm, Swamlnathan and Lewis (1984) 
and Elm and Lewis (1985) are some of these studies. Hlner 
and Lewis (1975) studied the body at angle of attack and 



restricted the analysis to the windward plane of 
symmetry. Kim and Lewis (1965) simulated the three 
dimensional flow field over the Shuttle orhiter. The 
analysis used shock and wall slip conditions with finite 
wall catalyticity. These works confirm the finding of Moss 
(1971) regarding non-equil ihrium heat flux predictions 
being lower than both the frozen and equilibritun flow heat 
flux predictions. 

Prabhu, Tannehill and Harvin (1987) solved the 

parabolised Havier-Stokes (PHS) equations for hypersonic 

laminar flow of a multi -component chemically reacting gas 

mixture over two-dimensional/axisymmetric bodies. Space 

% 

marching technique is used to evaluate the solution at the 
down stream stations. The equations are solved from body 
to the free-stream. The marching procedure requires that 
the local Mach number be greater than unity in the 
inviscid part of the flow. This implies that some other 
procedure is required to generate the solution in the 
stagnation region of the blunt body. The analysis shows 
higher degree of dissociation near the wall for PHS code 
in comparison to bOTindary layer results. This is due to 
higher pressure and density at the wall for PHS 
predictions. 

Bhutta and Lewis (1966) solved three-dimensional PHS 
equations for chemically reacting gas mixture over a 
sphere-cone body at large angle of attack. The analysis 



r Ecvi,-, £~ec:es £.riOL reaction nctce2. A nav 

v: vcicic-r -c-cr:‘fcclcr scnerje is devised to account for 
sirong cross-coupling in tne cross flow separatee regions 
on tne leeward side of tHe Icody. The how-shocK. shape and 
its location are computed as the part of the solution. The 
fluid-mechanical equations and species continuity 
equations are decoupled In the solution procedure. Their 
coupling is accounted through an Iterative procedure. 
However, this method can not he used to simulate the flow 
field in the stagnation region of the body. The analysis 
shows that with angle of attack, the leeward side has more 
dissociation of oxygen and nitrogen and formation of NO'*' 
than on the windward side. This is due to the thicker high 
temperature region on the leeward side. 

Park (1965) in his paper discusses the various 
aspects connected with the conqputatlon of the chemically 
reacting flow field. This work also discusses the 
Situation when vibrational temperature and translational 
temperature are not in equilibrium. The influence of this 
phenomenon on the reaction rates is also discussed. 

Park identifies the numerical problems encoiintered In 
chemical non-equilibrium flows as follows: 

i 

( I ) 'Violation of overall mass conservation, 

( II ) Instability caused by spiirious negative 
concentrations. 


Violation of mass conservation also leads to 
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of conscsnat. :>n of CMentum SiTxC. c-iiC-g: 




"Lc se pr-oMeris are fiiscus’sed by ParKL in connection witSi 
tbe unstea^Jy approacli, they are equally relevant to steady 
state solution procedures such as the viscous shock layer 
isethod. 

The two thin layer model predicts the stagnation point 
heat transfer rates in good agreement with the 
eKperlmental results and is in several respects superior 
to the higher order boundary layer theory. However, one 
should not expect the two thin layer theory to predict the 
detailed flow field correctly in the highly rarefied flow 
regime, due to several simplifying assumptions, which are 
not valid under rarefied conditions. These assumptions are 
as follows: 

(1) The basic assumption of zone of disturbance being 
thin and thus the order of magnitude analysis is not valid 
under rarefied conditions, where the shock wave and shock 
layer are thick in comparison to body nose radius. The 
terms dropped based on this ass\imption are not negligible. 

(il) The assumption of thin one dimensional shock is not 
valid and the tangential transport of mass, momenttun and 
energy across the shock can not be neglected. 

(ill) The coupling between the shock layer structure and 
the shock structure is very weak. The shock layer solution 
caul be obtained independert of the shock structure, in 
case of low density flows, there will be stronger coupling 


,c? f.r i-r^ocK struct'-re a::i the shock structvrer, 

a.^d the s2 qcK 2:trucf..:-«e will influence the shock layer 

iGiiruc vUre* 

The failure of the two thin layer theory has been 
attributed by Ahouse and Bogdonoff (1969) to the failure 
of Havier-Stokes e<Iuations in the rarefied regime. But 
this is not correct as this model does not represent the 
full Havier-Stokes equations. It is possible that 
considering the full Havier-Stokes equations in this model 
Will lead to better results and higher range of validity 

for this model. It is one of the aspect that is studied in 

detail in the present dissertation. 

Huroerical Integration through the Shock; 

In this approach the governing equations, either 
the full Havier-Stokes equations or a simplified form, are 
integrated from the body surface to the free stream 
through the shock. Depending on the flow parameters like 
Kach number and Reynolds number, the position of the sharp 
shock or the merged profile of a thick shock wave and a 

thick viscous layer near the body come out as an Integral 

part of the solution. However, the difficulties associated 
with the numerical integration of the Havier-Stokes 
equations make the computations very expensive. Also at 
higher Reynolds number, integration through the thin shock 
transition zone causes numerical instabilities. 

Levinsky and Yoshihara (i96£) and Eao( 1964b) used this 



approach in their investigation of the l>i'ant hody prohieis. 
Levinsky and Yoshihara obtained the exact numet'ical 
solution of sli^pllfled Havier-Stokes equations ' for 
Insulated and cold wall cases for a monatomic gas.^ In 
this investigation, thin layer assumption was used to 
simplify the Havier-Stokes equations, hut this analysis ■ 
contained more terms than originally used by Cheng in' his 
thin layer model. This extended the validity of the 
results to much lower Reynolds number than Cheng's theory 
in predicting the detailed flow structure. At Reoo = 152, 

it is noted from this analysis that the region of - 
disturbance is approximately 50Z of the body nose radius. 

Kao (1964b) used the full Havier-Stokes equations 
with local similarity assumption to investigate the 
stagnation region of a blunt body. Jain and Adimurthy 
(1974) extended this analysis further and obtained the 
solutions' for various free stream conditions. These • 
results show that the failure of the thin layer model to 
predict the correct behaviour of surface pressure -with' 
decreasing Reynolds number was not due to the failure of 
the Havier-Stokes equations, as envisaged by Ahouse and 
Bogdonoff (1969) and Vogenitz and Takata (1970) but was due 
to the simplifications introduced in the analysis by 
virtue of various assung>tions. 

Bellinger (1971) and Kumar and Jain (1975) 
obtained the merged layer solutions for cheiBlcaUy- 
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reacting flow in tlae stagnation region of a bi;mt ic /. 
However, both these analyses assume that viLveln. 
degree of freedom is fully excited and thus do not 
consider actual variation of the thermal energy of the gas 
as a function of the gas temperature. Also these 
investigations assume a simple power law variation of 
viscosity instead of a more realistic variation depending 
upon the collision cross-section. Surface slip effects are 
also ignored in the investigation. These studies show that 
the merged layer analysis predicts higher temperature in 

the shocK layer than the thin layer theory. The higher 

temperature gives rise to higher reaction rates and 
consequently higher levels of dissociation and ionisation. 
Gupta and Simmonds (19d6) extended this analysis by 
considering the actual variation of specific heat and 

viscosity of gases. This analysis employs the multi- 
component surface slip conditions. The energy equation Is 
used in terns of temperature, thus avoiding the use of 
time consuming process of converting the enthalpy into 
temperature for use in reaction rate estimation and 

transport property estimation. These results show a good 
agreement with Direct Simulation Konte Carlo (DSHC) results 
of Hoss and Bird(l965) for the case of shuttle orbiter. 

Several authors have used unsteady approach to solve 
the full Havier-stoKes equations for simulating the flow 
field over a given body. In this approach, the unsteady 
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equations are solved and time marcHing procedure Is used 
to get tne steady state solution. HacCormacK (1969) 
proposed liis predictor-corrector sclieme for solving tne 
full Havler-StoKes equations. The scheme is second order 
accurate in space and time. The constraint on the time 
step has also heen derived using stability considerations 
for the explicit scheme. However, due to very small time 
step, the time involved in the computation is very large. 

Widhopf and Victoria (1972) used unsteady approach 
for solving the full Havier-StoKes equations for 
chemically reacting gas mixture past a blunt body for 
merged layer flow analysis. The analysis used seven 
species model for the air. The comparison with thin layer 
solution shows that thin layer model predictions for 
concentrations of O, H and HO are lower as compared to the 
merged layer model using full Havler-StoRes equations. 
However, this difference between the predictions by two 
methods decreases with increase in Reynolds number. 

Li (1987) used unsteady approach to solve laminar 
chemically reacting flow over a blunt body under 
hypersonic low density conditions. The analysis used the 
full Havler-StoKes equations between the body and the 
shocK. A bow shocK is treated as a con^utational boundary. 
ShocK slip and surface slip conditions were not considered 
in the analysis. In the solution procedure, fluid 
mechanical equations and species continuity equations were 


decoupled. The approacla of decoupling tiie two sets of 
equations results in faster convergence. 

Candler and HacCormacK <1988) solved fhe full Kavier- 
StoKes equations for multi -coi^ponent gas mixture in 
tliermo -chemical non-e<iuilihri\am for a two-dimensional 
configuration. Implicit Gauss-Seldel line relaxation 
technique is used for solving unsteady Navier-StoKes 
equations to obtain steady state solution. The analysis 
uses three temperature model for energy conservation and 
estimation of chemical reaction rates. The analysis 
considers both five species and seven species model for 

t > 

the air. It is noted that stagnation heat transfer rates 
for reacting flow are nearly SOX higher than the 
corresponding perfect gas value. The analysis shows that 
electron temperature is approximately four times the free 
stream tenperature. Translational temperature is 
approximately twice the vibrational temperature except 
near the wall, where all the temperatures are nearly equal 
to the wall tenperature. 

Direct Slmalatiop Honte Carlo <PSHC} Technique 

There have been attempts to simulate the flow field 
using Boltzmann equations. In this approach, the fluid is 
not considered as a continuum. Particles are considered 
individually or in small groups while slmalating the flow. 
This approach is essentially valid for highly rarefied 
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flow regime, wlsere tlie Enudsen number is of oritr os 
iinity or liigner. Tne only closed equation tha-, ir. 
applicable is tbe Boltzmann equation. Solutions to tlte 
Boltzmann equation are readily obtained in tbe :.'c.e 
molecular limit as tbe Enudsen number tends to infinity. 
However, as tbe Knudsen number decreases and tbe 
transition flow is encountered, it becomes difficult to 
solve tbe Boltzmann equation. Honte Carlo procedures brve 
been developed to solve sucb problems. Bird (1966) used 
tbis method to evaluate beat flux and drag of a sphere in 
transition regime. Vogenltz and TaKata (1970) computed tbe 
flow-field in tbe stagnation region of a blunt body. Bird 
(1976, 1961) bas developed one sucb Honte Carlo technique 
called Direct Simulation Honte Carlo (DSHC) for simulating 
tbe transition regime flow. 

In tbe DSHC method, molecular motion and 
inter-molecular collisions are simulated using appropriate 
probability law. Tbe time parameter in tbe simulation can 
be identified with tbe physical time in tbe real flow, and 
all calculations are unsteady. For steady state flow 
conditions, time asymptotic limit of tbe unsteady soluticn 
is obtained. Tbe computation is always started from an 
initial state that permits exact classification sucb as a 
vacuiim or uniform equilibrium flow. Tbe computations laRt 
advantage of flow symmetries to reduce tbe number of 
position coordinates and tbe dimension of cell networX, 
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i:V-..o?r.^ For sirjjulatir-t: the rr. 're realistic viscous 
itehav-our of the fluid, variable Hard sphere (VH3) model is 
used. The V;1S model accounts for more realistic 
collisions between the molecules by modifying the hard 
sphere model. The chemical reactions 2 u:*e incorporated into 
the analysis by considering the reactive cross sections. 

Though the DSHC model is basically developed for the 

transition flow regime, currently it is being used to 

simulate the flow for the £nudsen number of the order of 

0.001. Moss and Bird (1985) used the DSHC technique to 

simulate the flow over Space Shuttle orbiter nose during 

re-entry for Knudsen number as low as O. 03. These results 

agree well the measured heat flux data in flight. Hoss, 

Bird and Dogra (1988) have applied the DSHC method to 

simulate the flow field during the atmospheric entry of 

Aero-assist Flight Experiment ( AFE ) vehicle. In this 

analysis, apart from chemical reactions, radiation is also 

considered. It is seen that at an altitude of 90 Km 

radiative heat transfer is nearly 1 percent of the 

convective heat flux, vdiile at an altitude of 80 km it is 

as high as 15 percent. The radiative heat flux is higher 

at lower altitude because of the higher overall 

» 

temperature in the shock layer. 

3 Surface Slip Condition 


In the low density flow regime, slip velocity and 



?0 

tewperature I * " ~ry c ..r: '~ion S Cj.X'ii j.' *. i£!, 'trf 

rarefaction ci-.'t >.£, Tvs vo -^h? low t,*:!!. j.ty ti_e tiov-', 
•(';.>. ■-' pc, x-ticies do noo sutfe.r- sufficient nureier of 

col 1 1 ^ n 4:is so as to lose tlieir entire momentum and tnus 
•,.of. c^cs a finite velocity relative to the solid surface. 
Siriiiarly, the temi erature of the gas adjacent to the body 
surface differs from that of the body. These boundary 
condit.ons become more important as the gas becomes more 
ar.d mere rarefied. 

Several authors have derived the expressions for slip 
velocity and temperature Juiqp boundary conditions using 
approximate solutions of the Boltzmann equation in the 
vicinity of the wall. These relations are derived under 
the assumption of small departure from local thermodynamic 
and chemical equilibrium conditions. The expressions given 
by ShidlovsKiy (1967) are; 
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viccr e j: and y are the coordinates .along and normal to the 
suric.ee respectively, while u and v are the tangentfal and 
r.orn;,: 1 velocities. Q and o are the momentum and thermal 
acc :r/„odation coefficients vf the surface and Ag is the 
ttisar. free path. It is to be noted here that even for 


coi^plete acconamodation (© = 0^1), ^.he slip *-cicr:,ty 
and temperature j'ump do not vanisli. Tne xc-rr- i:.vci\ jr.c 
temperature gradient in tlie expression for velocity olip 
is referred to as the thermal creep. According to this 
term, a temperature gradient along the surface induces a 
flow in the direction of increasing temperature, usually 
this term is neglected in the studies involving slip 
velocity as this terms belongs to a higher order of 
smallness than the velocity gradient term. 

In case of a multi- component gas mixture, the 
diffusion of species influences the slip velocity and 
temperature Jump conditions at the wall and thus needs to 
be accounted for in the distribution function. Scott 
(1973) evaluated the wall slip conditions for a multi- 
component gas mixture with diffusion and wall catalysed 
atom recombination. These Conditions were simplified by 
Hendricks (1974) for direct application to the fluid flow 
problem. However, these conditions led to an Enomalo'uS 
situation where temperature Jump and slip velocity 
decreased with increasing rarefaction ( Johnston and 
Hendricks 1978). Gupta, Scott and Hoss (1985) evaJuateci 
these conditions once again for multi -component gas 
mixture using a distribution function, which account'eo fcr 
the diffusion of the species more accurately. 1i:c£e 
conditions when applied to fluid flow pi^obleai oid not 


exhibit anomalous behaviour. Present analysis uses 


these 
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slip conditions whicli are reproduced in section 2.2. 

Tne derivation of the slip conditions assumes small 
departure from local thermodynamic and chemical 
equillbri'um conditions, that is, gradients of velocity, 
temperature and concentration are small. The validity of 
these slip conditions is not established for large 
gradients. 

1.4 Scope of the Present WorK 

In this worK, hypersonic flow in the stagnation zone 
of the blunt body is investigated for perfect and 
chemically reacting gases. In the previous sections 
several limitations of existing continuum approaches have 
been pointed out. Some of the unresolved questions in this 
regard are as follows: 

(i) Is it possible to extend the range of validity of the 
shock layer method by use of full Navier-Stokes equations? 
(li) How do the chemical reactions within the shock 
transition zone influence the shock layer flow-field at 
low Reynolds numbers, where the shock wave is thick? 

(Ill) Is it possible to extend the range of Reynolds number 
for which non- equilibrium computations converge? 

(iv) How to overcome the violation of mass conservation 
law normally occurring in numerical studies of chemically 
reacting flow? 

(V) Is it possible to overcome the numerical instability 



caused i>y spurious negative concentrations? 

Tne aim of tne present investigation x& to find oit 
the answers to questions posed above. In the present 
analysis full steady state Havler-StoKes equations with 
species continuity equations are assiuned to he valid for 
describing the hypersonic rarefied flow field past a blunt 
body. The region of interest is the stagnation region, 
which is most crucial for thermal protection system design 
considerations. In this region, the non-dimensional ised 
governing equations are reduced to a set of ordinary 
differential equations with split boundary conditions 
under the assumption of local similarity. Post-shocK 
transport effects are accounted for by using modified 
RanKine-Hugoniot conditions. Surface slip conditions for 
multi -component gas mixture are used at the wall. The 
analysis considers a real gas with variable specific heat. 
Transport . properties are evaluated using available 
correlations for multi -component gas mixture. The 
numerical solution to the governing equations adapted to 
the stagnation region is obtained using a finite 
difference procedure. 

As . a result of present studies, it is observed that 
the use of the full Havier-StoKes equations within the 
viscous shock layer approach provides uniformly accurate 
results both for high as well as low Reynolds niuuber. At 
high Reynolds number, heat transfer results approach the 


boundary layer values while at low Reynolds nuiid>er it 
gives correct behaviour of surface pressure and adiabatic 
wall temperature. It is seen that wall pressiu'e and 
adiabatic wall temperature increase with decreasing 
Reynolds number as observed experimentally. The present 
study shows that the use of full Havier-StoKes equati.ns 
within shock layer model increases the range of validity 
of the solution without having to solve fully merged layer 
equations from the body to the free stream. 

For the case of a real gas, the flow-field obtained 
from present computations agrees well with merged layer 
and DSHC results. It is also observed that in Addition to 
the Reynolds number, free stream velocity influences the 
heat flux to the wall for the real gas flow. It is also 
observed that chemical reactions within the shock 
transition zone do not appreciably influence the shock 
layer flow field. The proper treatment of the production 
terms in the species continuity equations has inproved the 
convergence characteristics of the solution procedure. It 
is also observed that non- equilibrium computations capture 
equilibrium concentrations in most of the shock layer 
regin at higher Reynolds number. 


The 

mathematical 

formulation of the 

problem 

is 

presented 

in 

Chapter 

Z. The solution 

procedure 

is 

discussed 

in 

Chapter 

3. In this chapter, 

details i 

about 


finite difference scheme, grid clustering and treatment of 
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mass production terms due to cliemlcal reactions are 
presented. At lilgher Reynolds number, grid clustering 
plays an Important role In resolving tbe details of tbe 
flow-field. Treatment of tne production terms also plays a 

crucial role In tlxe evolution of various profiles in the 

\ 

flow-field. 

Chapter 4 presents the detailed results obtained from 
the computations and their comparison with available 
theoretical and experimental results. The results are 
computed for a wide range of Reynolds number extending 
from the boundary layer regime on one side to the 
transition regime on the other. Detailed flow profiles as 
well as surface characteristics are presented. 

In Chapter 5, the Inportant conclusions drawn from 
the present study are' reiterated and summed up. Thus, this 
chapter focuses the contribution of the present study to 
hypersonic, chemically reacting, rarefied blunt body 
problem. A brief statement on the further research that 


can be carried out is also Included. 


CMAFTER 2 


MATHEMATICAL. EOEMUL AT I ON 

2. 1 IBTRODPCTIOH 

In the present study, the full steady Havler-StoKes 
equations for chemically reacting flows are employed for 
the analysis of the flow field In the stagnation zone of a 
blunt body. The equations are specialised for stagnation 
zone by expanding the flow variables about their 
stagnation line values using flow field symmetry In the 
stagnation zone of a blunt body. In order to account for 
the rarefaction effects present at low Reynolds ntunbers, 
surface slip conditions for multi -component gas are 
employed. Similarly at the shock, transport effects are 
considered while evaluating post- shock conditions. 
Section 2. 2 describes the governing equations and boundary 
conditions as applicable to the stagnation zone. 

For chemically reacting flows it Is necessary to 
consider various reactions occurring in the flow field and 
their Influence on the flow field through the mass 
production terms. The mass production terms depend on the 
reaction rate constats for various reactions. Section 2. 3 
describes the different reactions considered in the 
analysis of the flow field and their reaction rate 


a? 


constants. 

The analysis of real gas flow field re<xulres the 
estimation of thermodynamic and transport ' properties of 
the gas mixture as function of local composition pressure 
and temperature. Section 2. 4 describes the method used for 
estimating the thermodynamic and transport properties of 
individual gases, it also describes the method for 
evaluating the properties of the gas mixture from those of 
the Individual gases. 

g. g 60VERH1BG EQPATIOBS 

The equations employed in the present analysis 
for the simulation of flow field in the stagnation region 
of a blunt body are the steady Havier-StoKes equations for 
compressible reacting mixture. The governing equations for 
laminar flow of i multi-componemt gas In/ chemical non- 
eqillbrlum are available in various references namely, 
Hirschfelder, Curtiss and Bird(1954), Bird, Stewart and 
Lightfoot (i960) and Williams (1965). Recently these 
equations have been derived by Lee (1986) accounting for 
thermal non-equilibrimn between various temperatures. The 
various ass\imptions used in the present analysis are as 
foil ows : 

1. Flow is steady, laminar and axially symmetric, 
g. There is only one temperature, that is, there is 
thermal equilibrium. 

3.. Radiative heat transfer has been neglected. 
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4. Mass diffusion is binary and due to concentration 
gradient only. Hass diffusion due to teiqperature and 
pressure gradients is neglected. 

6. The influence of production rates for various species 
in tbe sbocK transition zone on tbe post<-sbocK conditions 
is neglected. Tbe post-sbocK transport effects are 
accounted for. 

7. Ambi -polar diffusion is considered for the charged 
species. 

The governing equations have been written in surface 
oriented coordinate system ( Fig 2.1 ). In case of flow 
around blunt bodies, the flow field is symmetric about the 
stagnation line. The flow variables have been non- 
dimensional ised as follows: 


u s u*/U«,* 
g : s / Rb 
P = p* / P** 

T s T* / T*j»ef 

US U* / T*pef) 

Re = p«* Rb/ u»( 


V = V*/U«»* 

= n / Rb 

P = P* /(P«* 2) 

h = h* / 

T*ref = / cp«, 

T*ref) 


( 2 . 2 . 1 ) 


The governing equations in terms of these non- 
dimensional variables can be expressed as follows, with ^ 
denoting the body slope and subscripts g and ^ denoting 
differentiation with respect to g and ■n respectively: 
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GlcDal ContlPuitT Equation 

( p ru)^ + tpr<l + ^ ) v).,,sO 

( 2. 2. a ) 

Tangential Hoaientian Equation 
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( 2. 2, 3 ) 
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Hormal Homeattam Bgaatlon 
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+ pw,^ - 
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(i+'n) 


p-n 
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( 2. 2.4 ) 


Energy Equation 
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P Cp ( 


Te + V T.,, ) - ( 


(I+'n) 
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Re r(i+‘n) 
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Evecins CoBtinmty Squat ion 


Cig + p V Ci^ = Wi 

1 +*n 

1 r M> Le r p. Le r(l-<-<n) 

+ ( i c ) + ( i ) 

Re rd+iD*- Pr (l+tj) ^ Pr "n 


for 1 = 1.2 H 5 

( £. 2.6 ) 


Egnatioa of State 


p T R 

p s 

CPoo H 

< 2. 2. 7 ) 

Here tue specific Heat Cp of tHe gas mixture is given Hy 


Hs 

Cp s E Ci CPi 
i=i 

( 2. 2. 6 ) 

and molecular weight H of the gas mixture is given by 


r Hs c - 1-1 

H = E i — 

«- 1 = 1 Mi J 

( 2.2.9 ) 

Bo tmdsry Conditions ; 

The governing equations ( 2. 2. 2 — 2 . 2. 6 ) along with 
supplementary equations ( 2 . 2,7 — 2.2.9 ) require 

boundary conditions wbicb are to be used while solving 
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tlaese equations. In the present case, as the solution is 
■ being sought between body and shocK, conditions need to be 
prescribed at the body and at the shock. Usually no slip 
conditions are applied at the body surface. However, in 
low Reynolds number regime, sufficient number of 
collisions do not occur so as to make the gas stick to the 
body and the flow slips at the body. The flow velocity and 
temperature of the gas at the body are not the same as 
those of the body. The surface slip conditions used in the 
present study for multi -component gas mixture have been 
derived by Qupta, Scott and Hoss (1965) and are given by: 
At the surface <n = 0 
Horiiial Velocity 

V s 0 • ( 2. 2. 10 ) 


Tanr^ential Velocity 
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( 2 . 2 . 11 } 
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Siirface TeiBPeratiire 
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Re p ^ 
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( 2 . 2 . 12 ) 


wJiere 


Ns Ns Ns H 

x= C Lei ^ t ^ J 

i=l j=l k=l Mu 


( 2. 2. 13 ) 


Ns H 1.5 

Y z E ( ) Ci 
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( 2. 2. 14 } 


Z = 1 + 1. 5 
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( 2. 2. 15 } 


Species Concentration 
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( 2. 2. 16 ) 
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witli 
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( 2. 2. 17 ) 
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< 2. 2. 16 ) 


However, tHough. equation (2.2.16) Is valid for 
general species concentration slip condition, for a non- 
catalytic wall one Has 


Cm * 0 

wblle for a fully catalytic wall 

Ci = Ci«, 


( 2. 2. 19 y 


( 2 . 2 . 20 ) 


THe slip conditions taken from Gupta et al. Have Heen 
adapted so as to avoid derived quantities sucH as y. 'IHese 
conditions reduce to no slip conditions as tHe Reynolds 
number increases to infinity. 
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stiecK Slip Ccndltlona 

The conditions at the shocK, nrhlch is being fitted in 
the present analysis, are normally given by KanKlne- 
Hugoniot relations. However, these relations are modified 
to account for transport effects down stream of the shocK 
at low Reynolds number. These shock slip conditions are 
given by 


p 'v r -Sin Ss 


u Sin Ss + — = Cos Sin Sg 

Re 


< 2 . £. 21 } 


( 2 . 2 . 22 ) 


p - V Sin Pg vfif s + Sin2 Pg 

3 Re p<jo» 


( 2. 2. 23 ) 


♦ V* 


Sin pg f h •(- 


li 

♦ Cp T4f = 

Re Pr 
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Sin pg { O. 5 ♦ ) + < 'll uifr + "v 'vA' ) 

Re 3 

Hs ij. Le^ 

+ E Ci^ hi 

i=l Re Pr 

( 2. 2. 24 ) 


and 



ti Le 






Ci * 

Re Pr 


( 2. 2. 25 ) 

Here ^ and’v are velocity components tangential and 
normal to tlie slxock and is tne outward normal to tHe 
snocK. Ps is tne shock slope. 

These shock slip conditions have been derived from 
Havier-Stokes equations under the assumption of thin one- 
dimensional shock. Under the assumption of thin shock, 
there is no contribution of mass production term to shock 
slip condition on the species continuity eqniatlon. These 
shock slip conditions reduce to usual conservation laws 
applied across the oblique shock with slope as the 

« Reynolds number increases indefinitely. 

Series Expansion 

The governing equations (2. 2. 2 - 2. 2. 9 ) and boundary 
conditions ( 2. 2. 10 - 2. 2. 25 } presented above can not be 
used directly to solve for the flow field at stagnation 
line. This is due to the term r, the distance from the 
axis of symmetry which is zero along the stagnation line, 
occurring in the equations. In order to specialise these 
equations along the stagnation line, the concept of local 
similarity is used and all the flow variables are expanded 
in terms of g about the stagnation line. It is further 
assigned that there is no influence of higher order terms 
on the behaviour of lower order terms. 
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Considering tHe synmetrr of the flow field about the 
stagnation line, the following series expansion is used: 


p s Po ♦ P2 . ( 2. 2. 26 ) 

p = PO + P2 , ( 2. 2. 27 ) 

T s To + Tg + ( 2. 2. 28 ) 

U S Ui g + Ug . ( 2. 2. 29 ) 

V s Vo { 1- g2/2 )+ Vg g2 + ( 2. 2. 30 ) 

Cj[ s ^10 * ^12 . ( 2. 2< 31 ) 

h : ho * hg I-. . , . ( 2. 2. 32 ) 

Pr = Pro + Prg g2 +. . . . <2. 2. 33 ) 

Le = Leo “*■ ^*2 6^ • ( 2. 2. 34 ) 

H = M-O + 1^2 • ( 2. 2. 35 ) 

Sin p .s Sin Ps S 1 - g£/2 ♦ ( 2. 2. 36 ) 

r s ( 1+ '*1) { g - g2/2 +. . . } ( 2. 2. 37 ) 

COS P = g - g2/2 < 2. 2. 38 ) 


Here the flow variables Po< Pgi Po> P2> P4 

I 

functions of <ri alone. It has been tacitly assumed here 
that the shock and the body are parallel to each other In 
the neighbourhood of the stagnation line. 

The series expansions ( 2. 2. 26 - 2. 2. 38 ) Is 
substituted in the governing equations ( 2. 2. 2 - 2. 2. 9 ) . 
Collecting the coefficients of like powers of g, one gets 
the governing equations which are applicable for the 
stagnation line. However, due to the elliptic nature of 
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tue Havier-StoKes equations, liigner order terms appear in 
tlie equations governing the lower order terms. All higher 
order temns except pg are dropped under the assiasption 
that higher order terms do not effect the behaviour of 
lower order terms appreciably. This assumption has been 
verified by E:ao(i964b) in the context of merged layer 
solution. Dropping the subscripts 0 and 1, one gets the 
following set of equations governing the stagnation line 
flow: 

Continuity Equation 

a p (u+v) +(1+11) ( p V = 0 

( a. a. 39 } 

Tangential Honentuni Equation 
p u (u+v) a ps 

+ p V u,^ + = 

1+1) l+'H 

ii 2 U« 8 U - 3v v<n 

1 + ^ - • 

Re 1+1) 9 (1+11)2 3(1+1)) 

mi u + v 

+ 1 I j 

ti 1+1) 


{ 2. 2. 40 ) 
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Horaial HoaentiHP Eguatiop 
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( 2. 2. 4i ) 


Equation for 


p u(u+v) 

Pa-n p V = 0 

1+11 


( 2. 2, 42 ) 


Energy Equation 
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< 2. 2.43 ) 
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fipccies Ccntlaulty eonatlon 
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i: 1. 2 H 5 (2. 2.44) 

The use of series expansions ( 2. 2. 26 - 2. 2. 36 } into 
the boundary conditions ( 2. 2. iO - 2 . 2. 25 ) gives the 

following boundary conditions: 


At the Surface 


•n = 0 


V = 0 


( 2 . 2.45 ) 
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Post- SHOCK Conditions 
At 'll s tje 

P ▼ s -1 


( 2.2.51 ) 


U U.,j s 1 ( 2. 2. 52 ) 
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( 2. 2. 54 ) 


M. Le 

Ci ♦ Ci^ = Ci« ( 2. 2. 65 ) 

Re Pr 


2 H 

P2 s -1 -V v<n ( 2. 2. 56 ) 

3 Re 


2. 3 CHBHICAL REACTIOB MODEL 

In tHe present analysis, a seven species and seven 
reaction otodel as proposed Hy Blottner (1969) Is 
considered. l!He species considered as constituents of air 
are Kg, Og, HO, O, H, HO^ and e". 'iSie model accounts for 
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a single 

Ionisation reaction 

only, namely, 

the 

Ionisation 

of HO. Eleven species 

model of 

Kang 

and 

I>\mn(1972) 

involving twenty six 

reaction 

was 

not 

considered 

in this study as It 

would involve 

more 

computation 

without significant 

Improvement 

in 

the 

confuted results as already noted by 

Swamlnathan, 

Elm 

and 


Lewis (1983). 

FIHITB RATE CHEMISTRY 

When the relaxation lengths associated with the 
chemical reactions are con^arahle with the flow 
dimensions, the rates of production of various species Wj^ 
are required to be evaluated. These terms appear in the 
energy conservation equation and the species continuity 
equation. For a multi -component gas mixture with Sg 
chemical species and H|{ chemical reactions, the chemical 
balance equation describing the overall change from 
reactants to products may be written as follows: 

HJ Kfy Hj 

£ ®lr ^1 <=ssssss> E Plr ^1 
Isl Kjjiv Isl 

N 

r = 1. 2 Hr 

< 2. 3. 1 ) 

Here Nj Is equal to the sum of the nun^er of the 
species and the catalytic bodies. The quantities and 

Pix* are the stoichiometric coefficients for the reactants 
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and tne products respectively for the reaction r, whereas 
^hr forward and backward reaction rate 

constants. Also, the quantities Xj^ denote the 
concentration of various species and catalytic bodies in 
moles per unit volume. 

The rate of change of concentration of any species as 
a result of given reaction r Is given by: 


dX 

( — 


at 


r 


( ^ 


ir 


-a 

ir 



HJ a HJ 

ir X - K IT 
j=l J j=i 


X^jr 

J 


• ( 2. 3. 2 ) 


This equation Is written in terms of X^, that is, the 
molar concentration per unit volume for various species. 
Writing It In terms of mass fractions of various 
species with the help of relations: 


Xi = p» Ci /Hi 

( 2. 3. 3 ) 

and 

dXi dp* 

Wi* = Hi ( ) = 

dt dt 

( 2. 3. 4 ) 

one gets 

Wip* = Mi { Pip - dir) 



Hj 

n 

J = 1 


c p* o 

< -1 ) 


Hj 

K V 
J = 1 






] 


i 2. 3. 5 ) 


6S 

This e<iuation gives the mass rate of production per 
unit volume of the species i due to the reaction r. The 
overall mass rate of production per unit volume for the 
species 1 is obtained by summing over all the 

reactions, that is: 

Hr 

Wi* = E Hi ( - Otij;») 

1=1 

HJ p» Hj T 

Z Tf ( — i - K ir ( — 1 — ) 
f*' J=1 Hj j = l Hj J 

( a. 3. 6) 

The seven chemical reactions considered in the present 
analysis are as follows: 

I 

r=l Og + Wi < > 2 O + u>j 

( 2. 3. 7 ) 

rs2 Hg + Wg < > 2 H + Wg 

( 2. 3. 6 ) 

r=3 Hg + H < > 2 H + H 

( 2, 3. 9 ) 

r=4 HO + <I >3 < > H + O + Wg 

( 2. 3. 10 ) 

r=5 HO O < > Og H 

( 2. 3. 11 ) , 

r=6 Hg + O < > HO + H 

( 2. 3, 12 ) 

. r=7 H + O < > H0+ ♦ e- 

{ 2. 3. 13 ) 

here w^, wg and Wg are the catalytic bodies. 
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Tbe reaction rates for these reactions are expressed in 
the modified Arhenlus form and are given hy 


Kfr = Cor T*Clr exp( - / T* ) 

( 2. 3. 14 ) 


and 

Khr * ^Or exp( - Dgr / T* ) 


{ 2. 3. 16 ) 


The constants occuring in these equations have been 
taKen from Blottner (1969) and are given in Table 1. This 
table also gives the third body efficiency for the gases 
considered in the analysis. 


2. 4 THBBMODTBAHIC AID TRAHSPOBT FKOFEKTIES OF GAS HIXTOBB 

Thermodynamic properties, specific heat Cpj, and enthalpy 
hj^^, are functions of teoqperatiLre for the species 
considered in the analysis. These are • required for 
estimation of post-shock conditions and estimation of heat 
flux due to diffusion of the species. Various transport 
properties like molecular viscosity , M-i*i thertnal 
conductivity Z*x and diffuslvity are also function? of 
ten^erattire and pressure. These are required for 

estimation of overall mixture viscosity, .thermal 
conductivity and diffuslvity. Since the multi -component 
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gas mixtures are considered to l>e mixtures of tliermally 
perfect gases, tne thermodynamic and transport properties 
for each species are computed using local ten^erature and 
pressure. The individual properties of these species are 
then used to evaluate overall properties of the gas 
mixture. Ihe following two sections describe the procedure 
for evaluating the Individual and mixture thermodynamic 
and transport properties. 

THERHOOYHAHIC FKOFEBTIBS 

The value of specific heat Cp^ and enthalpy h^ are 
obtained as function of temperature by using a fourth 
order polynomial curve-fit for specific heat Cp^ for each 
chemical species. The polynomial for enthalpy is derived 
from the expression for the specific heat. The polynomial 
for specific heat for species i is given by 

Cpi 

= a^ + ag T + a3 T*^ + ajj, T*’ + a5 t*^ 

R 

( a. 4. 1 ) 

and for enthalpy one has 

hj, ag T*2 aj a^ T*^ ag T*® 

s T* + + + * + ag 

R .2 3 4 5 

( 2.4. 2 ) 

where R is the universal gas constant. The values of 
the different constants of the polynomial (ai to ag ) are 
taKen . from ‘ Gupta (1957}. These constants have been 
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evaluated for each species for different tesqperature 
range. These expressions give specific heat and enthalpy 
per mole. The constants used In the present study for each 
species are gilven in Table S. 


TKABSFORT FROFEBTIBS 

The gas mixture is considered to be a mixture of 
thermally perfect gases. The transport properties of each 
species as a function of temperature are evaluated using 
Kinetic theory and later mixture properties are evaluated. 

, To a first approximation { first order Kinetic theory) 
the viscosity of a perfect gas( Hirschfelder, Curtiss and 
Bird 1964) Is given by 




2.6693x10-5 Hi ) 


<^1^ fill 


< 22 ) 


( 2. 4. 3 ) 

and the thermal conductivity for mono-atomic gases 
species is given by 


^*1 mono 


1. 9891x10“* 


< T^/Hi) 


0^2 Sii(£2) 


{ 2.4.4 ) 


or 




1 mono 

4 Hj, 


( 2. 4. 5 ) 


However, tne thermal conductivity for a polyatomic 
gas has an additional contribution due to the transfer of 
energy between the translational and internal degrees of 
freedom. This contribution to the thermal conductivity is 
given by 


2 CPi 

, 0. 88 < 1 ) mono 

5 Ri 

( 2. 4. 6 ) 

and thus 


16 Ui 2 CPi 

R ( 1 + 0. 88 ( 1)1 

4 Hi 5 Ri 


( 2. 4. 7 ) 


Here the different variables have the following units 
T* — temperature ( ^ K ) 

Hi — Holecular weight of the species i 

Oi — collision diameter of the species 1 < ) 

lii* — molecular viscosity of the species l (gm/cm sec) 

Zi* — thermal conductivity of species i (cal/sec 

fiii^S -- collision integral for viscosity 
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The collision integral is a function of the non- 

dimensional reduced temperature 




T+ 


(€l/K) 


( 4. 8 ) 

where is the maximum energy of attraction between 
the colliding molecules and K is the Boltzmann constant. 
The collision integral is approximated by VHhite 

(1974) as 


Qli£2s 1.147 { T+ ) -0.146 + ( T+ + 0. 6 ) “2 

( 2.4.9 ) 

The binary diffusion coefficient Dij is given by 


I>iJ 


2. 628x10"3 t* 5 

■ Hl*Hj 


2 Hi Hj 


p« aij2 fiijCli) 

{ 2. 4. 10 } 

where 

®lj^^ — collision integral for diffusion 
p* pressure in atmosphere 

’^Ij collision diameter for combination of species (A'^} 

Hi, Hj — molecular weights of species i and J 
I>lj diffusion coefficient { cm2 / ) 


ei 


Here for estimating the diffusion coefficient, the 
parameters and are evaluated as follows: 

dl + dj 
— 

( a. 4. 11 ) 

and 

eij = J<«i «J~)' 

( a. 4. la ) 

and the collision integral is approximated hy 

s ( + ( T* + 0. 5 )“2 

( a. 4. 13 ) 

with 

T+ = T» / <eij/k) 

< a. 4. 14 } 

The data regarding collision diameter dj^ and energy of 
attraction < e^ / k ) are given in Table 3 for the 
species considered in the present analysis. 


H1XT0BB FBOFBSTIES 


The mixture specific heat and enthalpy are given by 

Hs 

Cp = E Ci CPi 
isl 


( a. 4. 15 ) 


and 


Hs 

h* = E Ci hi* 
i=i 



< a.f. 1^ ) 


5e' 


vtiere C;|^ is tHe mass fraction for species i and Cp^ 
and Hi are specific Heat and entHalpy respectively for tHe 
species i. 

THe expression for viscosity and tHermai 
conductivity for tHe multl-con^onent gas mixture are given 
by HirscHf elder, Curtiss and Bird (1964) according to 
Kinetic tHeory. However, tHese expressions are ciuQbersome 
and computationaly time consuming. Therefore, most of the 
studies use an approximate technique for estimating the 
mixture viscosity and thermal conductivity. 

According to semi-en^irical formila of WilKe ( from 
BroKaw 1956, 1961), the mixture viscosity is given by: 


^*'*mix 


HS Xj[ 

B 

1=1 


Hs 


Where 


[ 




ij 


E 

6 

X 

i-i 

»ij 

3 



( 2.4.17 


0. 5 

H 0.25 -12 

1 + ( — 

) 

< 1- ) • 1 



Hi 3 

I — r 

M • 

,0. 5 

242 1 + - 

— i. 


L 

«j • 

j 


( a. 4. 18 ) 

Here ii“i» and Xi are viscosity, molecular weight and 
mole fraction respectively for species i. These relations, 
though developed empirically, have been ' Justified by 
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rigorous kinetic tneory ( BroKaw 1958). 

Tlie tnermal conductivity of a mixtiire of polyatomic 
gases can be divided into two portions 

j*^*mix ® ^*mix mono * *^'*mix 

( 2. 4. 19} 

here mono represents tne monatomic tnermal 

conductivity of tue mixture and gives tne 

correction due to transfer of energy between translational 
and internal degrees of freedom. 

Monatomic mixture conductivity is given by 


^*mix mono 


Ms 

E 

isl 


Xi 


Hs 



xj 


( 2. 4. 20 } 



where 


14 


r K* H n 

1 + ( L-jaans ) 0 . 5 ( _i ) 0 . 25 

^ i mttnet Hi 1 


[ 


4 mono 

M . <|2 

1 + — i- 

«4 


[ 


( H - H ) { H - . 142 M ) 
i ♦2.41 i i i 1 


c Hi + Hj )2 


< 2.4.21 ) 


and correction to the mixture conductivity 


K'Vx = 


Hs 

E 

isl 


g^*i »i 
Hs 

Six 
j=i J 


( a. 4. 22 ) 


where 




K* M n2 

1 mOQg ) 0. 5 ( _1 )0. 25 

mono 



( 2. 4. 23 ) 


However, in recent studies it has heen noticed that 
various assumptions regarding the interaction between 
molecules in the analysis of Wilke for mixture properties 
are not valid for temperatures above 10000 9 k. Armaly and 
Sutton (1980, 1981} considered these Interactions at these 
higher temperatures aind suggested an alternative 
relationship for mixture viscosity and thermal 
conductivity. These modifications in the correlations are 


as follows: 



is given by 



( 2. 4. 24 } 

where 

Aj^j = 1.1 for atom interacting with their own ion 
s 1,25 for all other interactions 
FiJ = i-0 

: O. 2 for atom and molecules interacting with 
electrons 

s O. 15 for atom and molecules interacting 
with ions 

= 0. 78 for atom and molecules interacting 
with atoms and molecules 
= 1.0 for all other interactions 

It IS worth noting here that this correlation of 
Arraaly and Sutton becomes identical to the correlation 
given by WilKe if one has =5/3 , = 1. 
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Tlie correlation parameter 


for mixture zht 


tij 

conductivity is given by 



1 

4^ Aij 


5 H 15 H 

+ 4 i- A + ( )2 

2 Mj 2 Hj 


( 1 + 


Hi 

) 

Hj 


2. 5 




JDQ^O^O. 5 
mono 



{ 2. 4. 25 ) 


Here the values of constants A^j , and F^j are 


K 


0. 25 B 


same 


as given in the viscosity model 


•taei_e: 1 


I : 


REACT I OH RATE CONSTAffTS 


Reac- 

tion 

Cor 

r 

Clr 

Cgr 

Dor 

Dip 

D£p 

1 

3. 61x10^® 

59400. 

-1. 0 

3. 01X10^5 

0 

-0. 5 

a 

1. 93X10^7 

113100, 

-0. 5 

1. 06X10^® 

0 

-0. 5 

3 

4. 15Xi022 

113100 

-1. 5 

a. 33X10®^ 

0 

-1. 5 

4 

3. 97x1020 

75600. 

-1. 5 

1. OlXlO^O 

0 

-1. 5 

5 

3. laxio^ 

19700. 

1. 0 

9. 63x10^^ 

3600 

0. 5 

6 

6. 75xlOl® 

37500. 

0 

1. 50xl0*i 

0 

0 

7 

9. 03x10® 

38400. 

0. 5 

1. 60X10*® 

0 

-1. 0 


THIRD BODY EFFICIENCY RELATIVE TO ARGON 


Efficiency relative to argon of 

Catalytic * 

h 'dies O ©2 H Hg HO 


Wi es 9 1 a 1 

v- 1 i - a. 6 i 

u'5 . ao 1 ao 1 ac 
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TABLE 2 

COHSTAHTS FOR POLYNOHIAL APPROXIMATION OF THERHODYHAHIC DATA OF GASES 


DATA FOR Ng 


Temp. 

ai 

aa 

as 


as 

ae 

300 

. 36748E0i 

i2081E-02 

, 23240E-05 

632i7E-09 

22577E-12 

1061 1E04 

iOOO 

. 28963E01 

. 15154E-02 

57235E-06 

. 99807E-iO 

65223E-i4 

905e6E03 

6000 

. 37270E0i 

. 46640E-03 

11400E-06 

. il540E-10 

32930E-15 

i0430E04 

15000 

. 96376E01 

25726E-02 

33019E-06 

14314E-10 

. 20332E-i5 

i0430E04 

25000 

51660E01 

23336E-02 

i2953E-06 

. 27872E-li 

2i359E-16 

10430E04 




DATA FOE 

: O2 


' 

Temp. 

ai 

aa 

as 

a4 

as as 

300 

. 36256E01 

i8762E-0a . 

70554E"-05 

67635E-08 

21556E-ii i0475E04 

1000 

. 36219E01 

.736ieE-03 

19652E-06 

. S6201E-10 

28945E-14 12019E04 

6000 

. 37210E01 

. 42540E-03 

28350E-07 

. 60500E-12 

5i860E-17 10440E04 

15000 

. 34866B01 

. 52384E-03 

39i23E-07 

. 10093E-il 

88714E-17 i0440E04 

25000 

. 3961 9E01 

. 39445E-03 

29506E-07 

. 73974E-12 

64209E-i7 i0440E04 
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TABLE 2 Continued 


DATA FOR HO 


Temp. 

ai 

ag 

as 

^4 

as 

as 

300 

.40459E0i -. 

34181E-02 . 

79819E-05 

-. 61139E-08 

. 15919E-11 

. 97453E04 

1000 

. 31890E01 . 

13382E-02 -. 

52899E-06 

. 95919E-10 

-. 64847E-14 

. 9828 3E04 

6000 

. 38450E01 . 

25210E-03 -. 

26580E-07 

. 21620E-11 

-. 63810E-16 

. 97640E04 

15000 

. 43308E01 -. 

58086E-04 . 

28059E-07 

-. 15694E-11 

. 24103E-16 

. 97640E04 

25000 

. 23507E01 . 

58643E-03 -. 

31316E-07 

. 60495E-12 

-. 40556E-17 

. 97640E04 



DATA FOE 

: O 




Temp. 

ai 

ag 

a3 

a4 

as 

as 

300 

. 29464E01 

i6381E-02 . 

24210E-05 

-. 16028E-08 

. 38907E-12 

. 29147E06 

1000 

. 25420B01 

27550E-04 -. 

31028E-08 

. 45610E-11 

-. 43680E-15 

. 29230E05 

6000 

. 25460E01 

59520E-04 . 

27010E-07 

-. 27980E-11 

. 93800B-16 

. 29160E05 

15000 

97871E-02 . 

12449E-02 -. 

16154E-06 

. 80379E-11 

-. 12624E-15 

.,E9iS0B06 

25000 

. 16428E02 

39313E-02 . 

29839E-06 

-. 81612E-li 

. 75004E-16 

. 29150E05 
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TABLE 2 Continued, 


DATA FOR H 


Temp. 

ai 

ag 

as 

a4 

as 

ae 

300 . 

25030E01 

21800E-04 

. 54205E-07 

56475E-10 

. 20999E-13 

. 56098E05 

1000 . 

24502E01 

. 10661E-03 

74653E-07 

. 18796E-10 

10259E-14 

. 56116E05 

6000 . 

27480E01 

39090E-03 

. 13380E-06 

11910E-10 

. 33690E-15 

, 56090E05 

15000 

12279E01 

. 19268E-02 

24370E-06 

. 12193E-10 

i9918E-15 

. 56090E05 

25000 . 

15520E02 

38857E-02 

, 32288E-06 

96052E-11 

. 95472E-16 

. 56090B05 




DATA FOR 

; H0+ 




Temp. 

ai 

ag 

as 

a4 

as 

as 

300 

. 36685E01 

il544E-02 . 

ei755E-05 

-. 48227E-09 

27848E-12 

. 11803E06 

1000 

. 28865E0i 

. 152i7E-02 

57531E-06 

. 10051E-09 

66044E-14 

. 11819E06 

6000 

. 22i41E01 

. i7760E-02 -. 

4303eE-06 

. 41737E-10 

-. i2828E-i4 

. ii619E06 

16000 

33240E01 

.24419E-02 -. 

19057E-06 

. 66580E-11 

-. 99112B-16 

. il819E06 

26000 

-. 43487E01 

.a4012E-02 

i4459E-06 

. 33813E-11 

28256B-16 

. il819E06 

4 , 


1 , ' 
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TABLE 2 Continued 


DATA FOR e" 


Temp. 

ai 

ag 

as 

a^ 

as 

ag 

300 

. 25000E0i 

0. 

0. 

0. 

0. 

74537E03 

1000 

, 25000E01 

0. 

0. 

0. 

0. 

74537E03 

6000 

. 25O80E0i 

63320E-05 

. 13640E-08 

10940E-12 

. 29340E-i7 

74500E03 

15000 

. 25001E01 

3ii28E-09 

. 35720E-i3 

i6036E-17 

, 25070E-22 

74500E03 

25000 

. 25001E01 

, 30157E-09 

22620E-13 

. 66734E-18 

66916E-23 

74500E03 





MOLECULAR CONSTANTS 

FOR TRANSPORT 

PROPERTIES 

>ecies 

e/K 

o 

H 


Ok 

ao 

g/g-mole 

O 

106. 7 

3. 050 

16. 00 

02 

106. 7 

3. 467 

32. 00 

H 

71., 4 

3. 296 

14. 00 

«£ 

71. 4 

3. 796 

28. 00 

NO 

116. 7 

3. 492 

30. 00 


NO-^ 


116. 7 


3. 492 


30. 00 


F|G:2.1 COOF4DS?4ATr SYSTEM 




OHAFTER 3 


METHOD OE SOLUTION 


3. 1 Introduction 

In order to get tne surface characteristics and the 


detailed 

flow 

structure, one 

has 

to solve 

the 

equations 

governing 

the 

flow 

field. This 

chapter 

describes the 

solution 

procedure 

adopted 

to 

obtain 

the 

numerical 

solution 

of 

the 

governing 

equations. 

The 

numerical 


solution involves the use of proper numerical scheme to 
ensure the convergence of solution. Moreover, in the 
present case the nonlinear production term becomes 
important and needs to he treated carefully for numerical 
convergence of the solution. Section 3. 2 describes the 
overall solution procedure, while sections 3.3 and 3.4 
describe the numerical scheme and treatment of the 
production term respectively. > 

3. 2 &0LUT10H FBOCEDUBE 

■The governing equations adapted for the stagnation 
ilo'vlf field were solved numerically by a finite difference 
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scheme. 

The set of governing equations contains Uq species 
continuity equations. However, three of these equations 
can he replaced hy algebraic relations using the condition 
of charge neutrality and elemental mass conservation. The 
condition of charge neutrality Is used to obtain the 
electron concentration and elemental mass conservation Is 
used to find the concentration of molecular nitrogen and 
oxygen. Thus only (Hg-S) species continuity equations need 
be solved. * , 

The set of equations Involves Hs second order 
equations namely, tangential momentum equation, normal 
momentiun equation, energy equation and (Hg-S) species 
continuity equations. Also It has 2 first order equations 
namely, global continuity equation and the equation 
for pg . There are ( 2 Hg+S) boundary conditions 
namely, tangential velocity, normal velocity, temperature ( 
or enthalpy ) and (Hg-3) species concentrations at the 
wall and also at the shock, and post shock values of 
pressure (or density ) and pg. Thus It Is seen that all 
the (2 Hg ♦ 2) conditions are required for solving the 
equations. 

However, in order to determine the sho6k stand-off 
distance one more condition is needed. As already noted by 
Ho and Probstein ii9W and Kao (1964b), this condition Is 
obtained by the application of global continuity equation 
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at the wall. At an impermeahle wall^the global continuity 
equation reduces to: 


This equation provides the extra condition needed for 
the detemination of the shocK stand-off distance 
However, the application of normal momentum equation at 
the wall results into: 


4 n r u 1 

p s — 21- <V -U) + V + 2V - — 2L_ -3. 5 u 

3 Re •- It 'H'H ^ 2 -• 

{ 3. 2. 2 ) 

If the normal momentum equation is solved for normal 
velocity V, it is imperative that the normal momentum 
equation he solved in such a manner that it satisfies not 
only two usual boundary conditions at the wall and at the 
shock but also the condition < 3.2.1 ) at the wall. At the 
same time, it is necessary that when continuity equation 
is solved for density ( or pressure through use of 
equation of state 2.2.7 ), it should satisfy the 

additional condition < 3.2.2 ) at the wall. Thus it 
appears that a second order differential equation is to be 
solved using three botindary conditions and a first order 
differential equation is to be solved using two boundary 
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conditions. 

In order to avoid tne difficulty of imposing tliese 
additional conditions, it was decided to solve the global 
continuity equation for normal velocity v and normal 
momentum equation for the pressure p. This avoids the 
apparent problem of imposing the additional constraints. 
Strictly speaking, these conditions have already been 
satisfied because these conditions were derived from the 
equations themselves. In this approach, the normal 
momentum equation, being first order in pressiire, requires 
only one condition and thus releases one condition ( post- 
shock normal velocity ) which can be used for determining 
the shock stand-off distance. An additional advantage of 

this scheme is that the solution approaches the boundary 

» 

layer solution when Reynolds ntuiber is very high. 

Once the free-stream conditions and the nose radius 
of the body are specified, initial profiles for all the 
dependent variables are assumed along with a guess value 
for shock stand-off distance. 

A unified finite difference representation of the 

« 

governing equations is generated as described in the next 
section, and the initial profiles are substituted into the 
finite difference form of the governing equations. The 
governing equations are solved in the order described 


below. 
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(i) Species Continuity Equations -- for mass 

concentrations of various species. 

(11) Energy Equation — for temperature . 

(ill) Equation of State — for density using old values of 
pressure and new values of temperature and mass 
concentrations of various species. 

(iv) Evaluation of tnermo dynamic and transport properties 
(V) Tangential Homentum Equation — for tangential 
velocity 

(Vi) Qlobal Continuity Equation -- for normal velocity and 
slxocK stand-off distance. 

(vll) Hormal Homentiam Equation — for pressure 
(viii) pg Equation — for pg 

The solution of the finite difference equations gives 
Improved profiles for the flow variables. This 

completes an iteration. These improved profiles are 
substituted into the governing equations once again 
yielding new profiles. This Iterative procedure is 
repeated till all the profiles converge to a prescribed 
accuracy. 

3. 3 Mumerlcal Scheme 

The tangential momentum equation, energy equation 
and species continuity equation are solved for tangential 
velocity, temperature and omss concentration of different 
species. All these equations are of second order with 



split boundary conditions. Usually in these eqpaations each 
derivative is replaced by corresponding finite difference 
representation. Thus these differential equations reduce 
to a set of difference equations which are to be solved 
for profiles of various flow variables. 

Choice of Grid 

In order to get realistic behaviour of the flow 
variables in the shocK layer, it is necessary that there 
are grid points at close intervals, wherever there are 
sharp gradients in the flow variables. If the grid is 
chosen to be uniform, large number of grid points will be 
needed to simulate the flow field and it will require 
large amount of computation time. In order to avoid large 
number of grid points and large amount of computation 
time, it is necessary that grid points are chosen to be at 
un-equal Intervals. In this way, one can have reasonable 
number of grid points ( and less coEQ>utational time) and at^ 
the same time get sufficiently accurate profiles for flow 
variables. 

In the present analysis, the grid chosen was 
geometrically stretched, with each succeeding grid spacing 
being 6 times the previous spacing as one proceeds from 
the body to the shocK. The first grid spacing is the 
function of 6 and nuifiber of grid intervals in the shocK 
layer. The grid stretching p-arameter <5 was chosen in such 
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a way tuat first grid point next to the wall was 
approximately 0. 5 percent of the boundary layer thickness 
away from the wall. 

However, if value of 6 is very large, then the grid 
spacing may be quite large near the shock. It is likely to 
lead to the breakdown of the solution procedure in the 
sense that the conditions at the shock may not have any 
influence on the solution near the wall. In order to avoid 
such a possibility, an upper limit for 6 has been 
prescribed for present analysis. This upper limit is a 
function of number of grid points enqployed in the 
analysis. 

Finite Plfference Scheaie 

For achieving smooth behaviour of the solution it is 
necessary that the system of difference equations be 
diagonally dominant. This condition is equivalent to local 
cell Reynolds number being less than 2 <Leonard 1979). 
This implies that the step size is quite small. This will 
result in large nuiid>er of steps and thus large computation 
time. 

One way to overcome the large number of steps is to 
introduce one-sided difference for- first order derivatives 
resulting in diagonally dominant difference equations. 
This ensures smooth behaviour of the solution variable, 
but it is at the cost of increased truncation error due to 
first order accuracy of one-sided difference 
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representation. In order to overcome this proDlem, in tne 
present analysis Unified Difference Representation (UDR) 
for differential terms was used. In this approacn DotH 
second and first order terms in a given equation are 
treated togetner and not as tne sum of two terms <Roaclx 
1982). Tne difference form is evaluated for botn tne terms 
comibined togetner. Witn the notations of tne Figure 3. 1, 
the difference form is given by: 


A fi^.ri + B = Cj fjL.j +C 2 fi +C 3 f^+i 


( 3. 3. 1 ) 


wbere 

Cl 


B [ i - exp<“BA2/A)3 

A 2 C exp(BAi/A) - IJ - Aj [ 1- exp(-BA 2 /A )3 


B C exp(BAi/A) -1 ) 

C3 = 

A 2 C exp(BAj/A) - 1] - Ai C 1- exp(~BA 2 /A)J 
Cg = -Ci -C3 


This form is always diagonally dominant for all 
values of A, B. A^ and Ag* For the case B= 0 , this expression 
3. 3. 1 reduces to the usual central difference form for 
second order derivative for unequal steps. Use of this 
scheme to second order differential equation results in a 
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set of tridiagonal difference equations can be 

represented for a grid point i as 

^i 2i-l + Zi + Ci ai+i 

( 3. 3. 2 ) 

wbere Zj[ is eitlier tangential velocity, temperature or 
mass concentration at tbe grid point i. gg^ 

tridiagonal system of difference equations was solved by 
forward and backward substitution scheme of Thomas 
algorithm. 

Qlobal continuity equation, normal momentum equation 
and P 2 equations are solved for normal velocity, pressure 
and P 2 respectively. These are first order initial value 
problems. These equations were integrated using 
trapezoidal rule. 

Continuity equation is also used to determine shock 
stand-off distance Hg. The continuity equation can be 
integrated to give 

1 2 fC 1 

P ▼ ( +C) = -2 ou({ + ) dC 

^e Jo ^g 

( 3. 3. 3 ) 

where C is the distance from the wall non-<li®®>^sio*^2Llised 
'aith respect to shock stand-off distance ^©* that is 



u.slng boundary condition at tne shock 
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at C s 1 p V s -1 


one has 

1 2 2 fi fl 

( 1 + ) - pudC+2 puCdt 

-He ‘He Jo Jo 

( 3. 3.4 ) 

This equation Is quadratic In and can he solved 
for- 'He' is noticed that for all the cases, one values 
of the nri^ Is positive and another negative. Thus there was 
no ambiguity In choosing 


BeMTKlaT»Y Conditions 

The governing equations are to he solved with the 
hotmdary conditions 2.2.45-2.2.50 at the wall. These 
conditions are mixed type, that is these involve both the 
variable and its derivative. In these equations, the 
derivatives are to he replaced by their finite difference 
foxvn and the resultant equations are used instead of 
governing equations at the wall. The condition 2.2.50 is 
used to evaluate the wall pressure from the pressure at 
the base of shock layer after solving the normal momentiun 
equation for pressure. Equations 2. 2.51-2.2.56 provide the 
boundary conditions at the shock. Ihese conditions are 
also mixed type and need replaceomnt of derivatives by 
their finite difference equivalent. 



However, wlien f ation is in terms of temperatiire, 
the post Shock condition given hy equation (2. 2. 64 ) in 
terms of enthalpy is in suitable fonn. It has to be 

rewritten in terms of temperature. The equation (2. 2. 54 ) 

can he rewritten as follows: 


ho + Cp ( T-Tq) + Cp T* 

Re Pr 


1 + 


1^* 


U<o 


»2 


v2 4 u Hs u Lei 

+ V v.r| ♦ E Ci* 

2 3 Re 1=1 Re Pr 


( 3. 3. 6 ) 

Here ho and To are the known post shock enthalpy and 
temperature at the previous iteration. T is the 
new post shock tenperature to he evaluated. This equation 
{ 3.3.5 ) is used in place of condition ( 2.2.54 ). 

3. 4 Treatment of Production Term 

It is observed from equations ( 2. 3. 6 ) that Wi* and 
hence Wi is a strong nonlinear function of Ci. This 
appears as a soxu'ce term in the species continuity 
equation emd in the energy equation formulated in terms of 
temperature. However, from the considerations of numerical 
convergence it is desirable to express Wi as a linear 
function of Ci. ihig can J>e achieved hy linearisation of 
term wi about the present value of Ci, Thus it is written 


in the fom 
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Wi = WiO - Wil Ci 

< 3. 4. i ) 

Similarly for energy equation also it is desirable to 
express tbe term E as a linear fvinction of 

temperatiire. It is written by expanding it about tbe 
current value of temperatiire in tbe following manner; 

Hs Hs k Ks dwx 

E bj^Wi = E ( biWi) + ( Tjj+i - Tjt) E bi ( )i^ 

i=i i=l isl dT 

( 3.4.2 ) 

Here k denotes tbe iteration at wbicb approximate 
profile is known and ( k+1 ) is tbe iteration for wbicb 
improved profile is being sought. Blottner <1969) bas 
discussed tbe need for splitting these terms in this form. 
Essentially this is required for getting tbe convergence 
of tbe numerical solution. If this expansion is not 
carried out tbeif one uses only the first, term of this 
expansion while evaluating tbe expression numerically and 
thus it amounts to neglecting tbe variation of this term 
with respect to or T during a particular iteration. 

If sufficient care is not excercised in tbe 
computation of chemically reacting flow, tbe net mass 
productiofa rate of all tbe species may not vanish (Park 
1985). Apart from the numerical errors, it also occurs 



because of tbe linearisation of tbe production terms < 
Appendix i). In order to avoid tbis problem, mass 
conservation for elemental species was used to evaluate 
tbe mass concentration of oxygen and nitrogen ( Coa and 
C|i2). It was also used to evaluate tbe production rate of 
oxygen and nitrogen ( woa and WHg) and tbeir temperature 
derivatives as well. 

For elemental mass conservation one bas 


Ho 

Coa t ^HO ^HO+ ) = 

' %0 


( 3 . 4 . 3 ) 


and 


Hu 

Ch2 ♦ Ch + ( Cuo ♦ Cho+ ) * Vu 

Huo 

^ < 3 . 4.4 ) 

Since production rate w^ is related to mass concentration 
by 

idCi 

Wi = p 

dt 


one gets 


Ho 

Woa * - Wo < Who ♦ who+ ) 

Huo 
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and 

Hh 

^H2 = - < Who + WHO+ ) 

Hho 

( 3. 4. 6 ) 

in the present work equations ( 3. 4. 5 ) and ( 3. 4. 6 ) 
are used to evaluate the production rate wog and wng 
use in equation ( 3. 4. 2 ) for energy equation solution. 
This ensures that at every Iteration net mass production 
for all species is zero and thus total mass is conserved. 

In order to evaluate terms dw^g/dT and dwng/^T for 
use in equation ( 3. 4. 2 ) , differentiating equations 

(3.4.5) and ( 3.4.6 ) one gets, 

dWo2 dWQ Hq 

dT . dT Mho 

and 

dWH2 dWH . Hh ewHo ^Who+ 

. ^ + j 

dT dT Mho dT dT 

( 3.4. 8 } 

Use of the equations ( 3. 4. 7 ) and ( 3. 4. 8 ) ens'ures 
that rate of change of overall mass Is zero with 
temperature. 

One more aspect of chemically reacting flows which 
' needs attention is the oscillations of profiles from one 


dWHo dWHo+ 

( + ) 

dT dT 

( 3.4.7 ) 



> tertion to another, '^is is possibly due to conversion of 
t:aergy from translational mode to chemical mode and vice 
ersa ( ParK 1966). Some times it may lead to one 
particular concentration becoming negative and this will 
further amplify and lead to other concentrations becoming 
i^egatlve, if unchecKed. Finally, it will lead to breakdown 
of the solution procedure. It has been taken care of in 
>be following manner: 

It is assumed that at any given point in the flow 
-i ield, the various ‘concentrations must lie between free- 
rtream value and maximum permissible values. 

Ciflo i Cl < Cl max 

for the species which do not have significant 
concentrations ih the free stream, namely O, H, HO, HO'*' 
and e~, usually referred to as minor species. The maximum 
permissible values can be evaluated from free stream 
conditions. In case concentrations during iterative 
procedtire are outside these bounds, then they are forced 
to the nearest bound. This ensures that none of the 
concentrations become negative during iterative solution 
procedure and eliminates problem associated with negative 
concentrations and unrealistically high values of 
concentration. Elemental mass conservation procedure is 
applied after all the above minor species have been 


evaluated. 



3> 6 Convergence of lumerlcal Procedure 


a ■ 


The procedure described above was coded in FORTRAH IV 
and was run on CDC Cyber 170/730 at VSSC Trivandman for 
evaluating the flow field and surface characteristics for 
various free stream conditions. The parameters used in 
connection with n\imerical solution procedure are as 
followst 

Humber of intervals : 100 
Grid stretching parameter : i.Oi < 6 < 1.05 
Convergence criterion : iX 

The convergence criterion of 1 percent means that all 
the variables over t^e entire shocK layer are within 1 
percent of their respective values in previous iteration. 

The computations were performed both for perfect gas 
and for chemically reacting non- equilibrium flow. For the 
perfect gas case, the free stream and fluid property 
parameters are given by: 

Re s 100 
M® = 10 
Pr = 0. 7 
y s 1. 4 
U » f T 

Figures 3. 2 - 3. 3 show the convergence of various 

profiles for the perfect gas case. Figure 3.2 shows the 
maximum fractional change for a flow variable over the 
entire profile. It is observed that maximum fractional 
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cbdiice In ' Ut v and T decreases to a value less . tban 
conveargence criteria by fifth iteration. 

Figure 3. 3 shows the evolution of temperature profile 
as l-terations progress. It is noted that at iteration 2, 
there is an increase in shocK stand-off distance which 
decreases as the profiles evolve with the iterations. It 
is observed that the profiles are more or less identical 
at it.eration 3 and iteration 4. Similar behaviour is noted 
in Figure 3. 4 which shows the evolution of tangential 
velocity profile. 

For chemically reacting flow, the free stream 
condlt-ions used correspond to the Shuttle flight 
conditions. These are as follows: 

Altitude s 92.35 Km 

) 

Free Stream Velocity s 7. 6 Km/s 
Hose radius Rb = 1. 296 m 
Wall Temperature = 1043 
Wall Catalyticity : Hon-catalytlc 

The thermodynamic and transport properties of the air 
were evaluated by the method described in section 2. 4. For 
mixture properties estimation, correlation proposed by 
Armaly and Sutton (1960, 1981) has been used. 

Figures 3. 5-3. 7 shows the maximum fractional changes 
in concentration profiles, flow variables and transport 
property with Iteration co\mt. The solution based on 
convergence criterion was obtained at twelfth iteration. 
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Fig\are 3. 5 shows the fractional change in the 
concentration profiles with iteration count. It is 
observed that masiiauin changes occur in concentration of 
atomic oxygen and therefore atomic oxygen profile is the 
last one to converge. This figure does not show the 
changes in electron concentration, molecular oxygen and 
nitrogen. The fractional changes in electron concentration 
is identical to fractional changes in concentration. 
The changes in concentration of molecular oxygen and 
nitrogen are related to‘ the changes in other 
concentrations by elemental mass conservation. 

Figtire 3. 6 shows the convergence of flow vaui'iables. 
It is observed from this figure that pressure is the first 
variable to converge, while tangential velocity gradient 
IS the last one to converge, i’igure 3. 7 shows the 
convergence of transport properties. It is seen that 
Prandtl number converges to the prescribed accuracy at 
iteration 2, while viscosity converges at iteration a. 

Figure 3. 8 shows the evolution of temperatxire profile 
with iteration count. It l^. seen that at iteration 1, 
there is a sudden decrease in shock stand-off distance. 
This is due to fact that initial guess value of shock 
stand-off distance correspond to the perfect gas 
conditions. It is seen that both post shock temperature 
and gas temperature at the wall evolve with the iterations 
and attain the final converged values. 
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Figures 3. 9-3. 12 show the convergence of 
concentration profiles with iteration count. Figru'e 3.9 
shows the evolution of concentration of molecular 
nitrogen. It is observed that with each iteration, level 
of concentrations decrease continuously and convergence is 
achieved without oscillations in profiles. Similar 
behaviour is noticed in figure 3. 10, which shows the 
evolution of concentration of molecular oxygen. It is seen 
that near the wall molecular oxygen is almost fully 
dissociated at about 2000 This is because near the 
wall, diffusion of species dominates over the chemical 
Kinetics. Figure 3. 11 and 3. 12 show the evolution of 
concentration profiles for atomic oxygen and nitrogen 
respectively. 

In order to assess the convergence characteristics of 
the present formulation and the computer code at higher 
Reynolds number, the flow field was computed for 
conditions corresponding to atmospheric entry of Aero- 
Assist Flight Experiment (AFE) vehicle. These conditions 
are as follows: 

Altitude s te. 90 Km 
Free Stream Velocity = 9. 365 Em/s 
Hose radius Rb s 2. 530 m 
Wall Temperature = 1111 
Wall Catalyticity : Hon-catalytic 
The thermodynamic and transport properties of the air 
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were evaluated, by tlie method described in section 2. 4. For 
mixture properties estimationi correlation proposed by 
Armaly and Sutton (1980, 1981) has been used. 

In this case the post-shock Reynolds number Re 3 is 
1760, while Hoss (1971) reported the failure of 
convergence of numerical computations for post-shock 
Reynolds number greater than 500. Moreover, in order to 
check whether the numerical solution has any tendency for 
divergence at later stages of the computations, the 
computations were carried upto 100 iterations, imen 
solution satisfied the convergence criterion at about 48 
iterations. 

Figures 3. 13 to 3. 18 show the convergence of various 
flow variables with iteration count. It is observed that 
xhe maximum change in all the variables over the entire 
shock layer decreases continuously after first 10 
iterations. Figure 3. 13 shows the maximum fractional 
change for both tangential and normal velocities and 
temperature. It is observed that these variables converged 
to the prescribed accuracy by about 30 iterations and 
there is monotonic decrease in the fractional change 
thereafter. Figure 3. 14 shows the maximum fractional 
change in the pressure and fractional change in the heat 
flux to the wall. It is noticed that pressure profile 
converges to the prescribed accuracy by about 30 
Iterations udille heat flux values satisfy this test after 
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al>out 10 iterations. This accuracy tei t has not hecn 
prescribed for heat flux value, nevertheless the heat flux 
to the body is found to have converged very fast. The 
convergence of heat flux by about 10 iterations implies 
that after 10 iterations, there are no significant changes 
in the flow variables near the wall. 

Figure 3. 15 shows the maximum fractional change in 
the concentration of Hg, Og and NO. It is observed that 
after initial sudden changes for about 10 iterations, the 
changes in the concentrations decrease monotonical ly with 
iterations. It is seen that concentration of Kg converges 
last at about 50 iterations. Figure 3. 16 shows the 
maximum fraction^tl changes in the concentration of O, H 
and HO'*’. The changes in the concentration of e“ is 
Identical to that of HO'*’. It is noted that in this case 
also fractional changes decrease continuously with 

iterations and concentration of NO'*’ is last to converge. 
Tne fractional changes in the flow variable do not show 
any tendency of increafi.ng with iteration except in the 
begirding when the flow variables are adjusting to -the 
boundary conditions*. Figure 3. 17 shows the convergence of 
shocK stand-off distance with iterations. It is observed 
that shock stand- off distance converges by about 25 
iterations and then it remains steady. Figure 3. 16 shows 
the convergence of of the wall heat flux. Heat flux 
converges by about 10 iterations and then remains 



ccr.stant. 




Figures 3.19 to 3. E9 sliow tHe evolution of the 
various p.^ofiles with n\imher of iterations. Figures 3. 19 
and 3. 20 show the evolution of the tangential and normal 
velocixy profiles. It is seen that velocity profiles 
converge hy about 30 Iterations and after 30 iterations, 
there are no appreciable changes in the profiles. It is 
observed that profiles at different iterations terminate 
at different distance from the wall as shocK stand-off 
distance evolves with the iterations. The profiles after 
60 iterations have not been plotted. 

Figures 3. 21 - 3. 23 show the evolution of 
temperature, c^nsity and pressure profiles respectively. 
In this case also it is noted that the profiles converge 


by 30 iterations. It is seen 

that 

temperature 

first 

increases 

with distance from 

the 

wall 

and 

later 

^€^Apel•ature 

remains constant. 

Near 

the 

shock, 

the 


temperature suddenly increases. This sudden increase in 
the temperature is due to post-shocK bounaary condition. 

Figures 3. 24-3. 29 show the' evolution of various 
concentration profiles. It is teen that profiles evolve to 
the converged profiles without oscillations. Figure 3. £4 
shows the evolution of Ng concentration K'oiile, It is 
observed that near the wall, initially concentration 
decreases due to dissociation of Hg as temperaiure 
increases. After some distance,- the concentration remains 



constant and men near me sliorK, me concentration of Hg 
increases sJiarpiy. This sharp increase in the 
concentration is due the boundary conditions at the shock. 
Similar behaviour is noticed in Figure 3.25, which shows 
the evolution of Og concentration profile. It is seen that 
at the wall Og is completely dissociated, while Hg is not 
fully dissociated. This is due to the wall temperature 
being low for Hg dissociation to take place. 

Figure 3. 26 shows evolution of the concentration 
profile of the HO. In this case it is observed mat near 
the wall HO concentration is nearly zero and the 

concentration increases with the distance from the wall as 
the temperature increases. Hear the shock, there is sudden 
decrease in HO concentration due me boundary condition. 
Figure 3. 27 shows the evolution of the concentration 
profile of O with iterations. It is observed that near the 
vail, concentration of O is large indicating complete 
dissociation of Og. This is due to fact that near the wall 
concentration is dominated by diffusion and not by 
chemical kinetics. Near the shock, concentration of O 
decrc ases suddenly due to the boundary condition. Similar 
phenomenon is observed in Figures 3. 28 and 3. 29, which 
show me evolution of concentration profiles of H and HO+ 
respectively. 

The computations were carried out for higher niimber 
cf intervals also, namely 200 and 600 for a few sample 
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cases. It was observed that the final converged profiles 
were identical n. th the converged profiles obtained for 
100 intervals. In view of the computation time involved, 
rest of the computation were performed for 100 intervals 


only. 
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RESULTS AND DISCUSSION 

4. i Introduction 

Based on tlie formulation given in cauapter 2 and tlie 
numerical procedure described in Chapter 3, a ntimber of 
computations were performed both for perfect gas and real 
gas for various free-stream and body conditions. These are 
discussed in the present Chapter. The basic purpose of 
these studies was: 

(i) To verify the conjecture regarding the failure of 
shocK layer approach, that is, failure is due to 
neglecting the terms in the governing equations and not 
due to the ShocK layer approach, 

(ii) To compare and bring out the differences between 
results from full Havier-StoKes equation solution and 
ShocK Layer model ( Davis 1970 ) with truncated Havier- 
StoKes equations and 

(iii) To address to the problem of numerical problem of 
numerical convergence of chemically reacting flow solution 
procedure including the violation of overall mass 


conservation. 



laa 

llie i^ftsults are l>elne presented in two parts. Section , 
4.2 discusses tne results for me perfect gas case wtille 
section 4. 3 presents results for me real gas case. 

4. 2 icsulti for Perfect flas 

Ttxe coBsmtatlons were performed for bom me Shock 
Laifer model < Pawls 1970 ) and for me present model of 
solving me full Havier>Stokes equation wimin shock 
layer, me free stream conditions were as follows: 

=.10 

Tw / TO s 0. £ 

Pr =0.7 

y =1.4 

U S f T 

e 3 Os 1.0 

The Reynolds nuimer based on reference conditions was 
varied between 100 and 1. These parameter were' chosen so 
that it is easier to compare me results wim reported 
results of Davis (1970). 

Comparison of overall caiaracteri sties 

The conflations were performed for me freestream 
conditions described above. For mese conditions it is 
noticed that me two approaches, namely, me present 
approach of using full lavier-Stokes equations and me 
Davis -approach of using truncated lavier-Stoluis equations, 
differ bom in overall characteristics and detailed flow 


/ 



field. TStese are discussed below. 

Figure 4. 1 presents the variation of shock stand-off 
distance as a function of Reynolds ntudber for the two 
methods. It is seen that the behaviour of shock stand-off 
distance is similar in the two methods. It is noted that 
initial 17 With decrease in the Reynolds number, the shock 
stand-off distance is nearly constant, but with further 
decrease in the Reynolds number, the shock stand-off 
distance increases. In the beginning, at higher Reynolds 
number, the viscous effects are confined to the region 
near the wall and do not influences the shock stand-off 
distance. Later with the decrease in Reynolds number, 
viscous effects encompass the* entire flow field and result 
in an increase of shock stand-off distance. It is to be 
noted here that the shock stand-off distance predictions 
by full llavl er-S tokes equations are lower than the 
corresponding predictions of the Viscous Shock Layer (VSL) 
model of Davis (1970). It was expected that the prediction 
of shock stand-off distance from the full Savier- Stokes 
equations will be higher than the VSL la'ediotions. The 
very assumption that the viscous effects are dominant in 
the region near the body ( argument used in order of 
magnitude analysis) implies that actually viscous effects 
encompass larger region though the viscous influence amy 
be small in this extra region, fhe higher prediction of 
shock stand- off distance by the vsL model is partly the 
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result Of tbe shock slip condition as applied in the VSL 
model. The shock slip conditions as applied in the VSL 
model are inconsistent with the govejmlng equations of the 
VSL model, under the shock slip conditions used by Pavls, 
post shock pressure is same as predicted by kanklne> 
Httgonlot relations. With shock slip conditions, the post 
shock temperature is lover than the temperature predicted 
by Ranklne-Hugonlot relations, which results in higher 
post shock density and lower post shock velocity, ^is 
lower post shock velocity will result in lower shock 
stand-off distance. Lover post shock velocity will in 
actual case result in higher post shock pressure and thus 
higher post shock density. This will result in further 
decrease in the shock stand-off distance, nie VSL approach 
of Davis (1970) does not account for this change in the 
post-shock pressure due to the shock slip conditions and 
thus predicts higher shock stand-off distance. 

Figure 4. 2 gives the variation of shock stand-off 
distance for the case of hot wall ( TwsTO ) with free 
stream conditions remaining the same, it is noted here 
that there is a continuous increase in shock stand-off 
distance with decrease in Reynolds number. In this case 
due to higher wall tengperature, the viscosity is higher 
and viscous effects influence a larger region of the flow 
field even at higher Reynolds number. 
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Flfure 4. 3 and 4. 4 sliov the variation of post shock 
pressure and post shock density respectively with the 
Reynolds nunher. It is noted that the present coonputatlons 
predict higher values of post shock lu^essure and post 
shock density as compared to the VSL computations. This 
behaviour is the outcome of the use of incorrect shock 
slip condition in case of the VSL solution as explained 
earlier. 

Figure 4. 5 shows the variation of pressure at the 
base of the shock layer, that is pressure of gas at the 
surface of the body, with the Reynolds number. It is 
noticed that at higher Reynolds number, present 
computations predict nearly constant base i^essure with 
decreasing Reynolds number, but with further decrease in 
the Reynolds number base jo^essure increases. This 
phenomenon is not observed in case of the VSL solutions 
where the base pressure continuously decreases with 
decrease in Reynolds number. However, wall pressure 
Initially decreases with Reynolds number and with further 
decrease In Reynolds number, the wall poressure increases < 
Figure 4.4). Vail pressure is computed from the base 
pressure using the shock slip condition for the px^sstu^e. 
This phenomenon of increasing wall and base pressure with 
decreasing Reynolds number is noticed in ease of the 
Havier-Stokes Solution only and not in case of the VSL 
solution. However, from experimental results ( Potter and 



Bailey 1964, Sherman 1963) it is Known that wall iwressure 
does initially decrease with decreasing Reynolds number 
and later increases with further decrease in Reynolds 
number. This clearly shows that full Sawier-Stokes 
equation solution procedure is superior to the VSL model 
in the low Reynolds number regime. 

It was stated earlier that the shock slip conditions 
used in case of the VSL solution are inconsistent with the 
governing equations and US's of proper shock slip 
conditions will lead to increased post shock ju^essure and 
thus overall higher pressure in the shock layer. In order 
to determine Whether the continuous decrease in wall 
pressure is due to inconsistent shock slip condition or it 
is inherent in the governing equations, cootputatlons were 
performed for VSL model with shock slip conditions 
consistent with the governing equations of the VSL model. 
Figure 4; 7 shows the variation of wall pressure with 
Reynolds ntimber as obtained from solution of Havier-Stokes 
equations and VSL equations with consistent shock slip 
conditions. This figure also shows that wall pressure 
continuously decreases with decreasing Reynolds number 
for VSL solution even with correct shock slip conditions. 
This shows that correct behaviour of wall pressure can not 
be predicted within the framework of VSL model and it is 
necessary to make use of full Havier-Stokes equations 
within the shock layer model to get the correct behaviour 
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of wall pressure wltli Reynol<^* Ji'Wober, 

Figure 4. 8 and 4. 9 si:iow tHe variation of lieat 
transfer coefficient wttn Reynolds numbej. for cold wall 
case ( Tw = 0.2 TO) and Hot wall case ( Tw = TO) 

respectively. For tne case of col<i case usual 

definition of tHe heat transfer coefficient based on 
tea«>erature differential Has been used while for the case 
of Hot wall Heat flux i* nondimensionalised with respect 
to Kinetic energy of tHe free stream as the temperatwe 
differential is zero. Figure 4. 8 sHows that the Havier- 
StoKes equations predict Higher Heat flux as compared to 
the VSL model. This is basically due to smaller shock 
stand-off distance leading to higher temperature gradient 
at the wall and difference in the temperature profile 
Itself. This figure also shows the results of Fay and 

Riddell (1958) and Lees (1956). It is seen that for Higher 
Reynolds number all the methods give results close to one 
another except ' Lees" method which does not account for 

variation of transport properties of the gas due to change 
in the wall teoaperature. With reduction in Reynolds 

number, both the numerical methods namely the present 

method of solving full Havier-stoKes equations and the VSL 
model predict higher heat flux than Fay and Riddell's 
method. This increase in heat flux can be attributed to 
the influence of the post -shock ^©rticity. Later on with 
further decreaste ®*y*ioids nuinher, surface slip 
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doou.na'tfts and Xiaat flux pradietlons toj numarical metlxods 
are lower la comparison to Pay and Riddell <1936) 
predictions. Plfiire 4. 9 snows tne variation of neat 
transfer coefficient wltn Reynolds number for tne ease of 
not wall. It Is noticed tnat Havler-StoKes equations 
predict positive neat flux to tne wall tnat is neatlng of 
tne wall for lower Reynolds number wnile VSL model 
predicts negative neat flux tnat is cooling of tne wall 
for tne entire range of Reynolds number, mis implies tnat 
according to tne solution of tne Havier-StoKes equations 
tne recovery temperature is greater tnan tne stagnation 
temperature of tne free stream. In case of VSL model, as 
neat flux Is always negative, it implies tnat recovery 
temperature is always less tnan stagnation ten^erature. It 
nas been noted from experimental data ( HlcKman and Gledt 
1963) tnat at low Reynolds numbers, tne z^covery factor 
increases beyond unity, mis mows tnat in order to get 
tne correct benavlour of beat flux wltn Reynolds niunber. 
It is necessary to consider tne full Havler-StoKes 
equations Instead of tne truncated Havler-StoKes 
equations. 

Comparison £f Structure Qi Plowfleld 

plgures 4. 10 tnrougn 4. 13 snow tne profiles of tne 
tangential velocity, temperature, pressure and density 
respectively for tne low Reynolds number case. 


fhese 
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figures also give the corresponding profiles for VSL 
model. It Is seen from Figures 4. 10 and 4. 11 that the 
variation of tangential velocity and temperature across 
the shock layer are similar in nature, hut gradient at the 
wall is higher in the case of Havler-Stokes equations 
solution. !thls results in higher skin friction and higher 
value of heat flux. 

Figure 4. 12 shows the variation of pressure across 
the shock layer. It is seen that apart irom quantitative 
difference in pressure, there is a qualitative difference 
in the variation of pressure in the shock layer. Havier- 
Stokes eqpiatlons not only predict higher post-shock 
pressure but also show continuous increase in pressure as 
one moves closer to the wall imile in ease of the VSL 
model, pressure is almost constant across the shock layer. 
This feature of pressure being almost constant across the 
shock layer is essentially a feature of the boundary layer 
and results from the fact that in the VSL equations all 
the viscous terms have been dropped from the normal 
momentum equation. Figure 4. 13 shows the variation of 
density across the shock layer. As expected, the behaviour 
of the density is different in the two models in view of 
the difference in the pressure profiles. 

Figures 4. 14 through 4. 17 show the profiles of 
tangential velocity , temperature, la'essure and density 
respectively at Reynolds number of lOO. it is seen from 
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Figures. 4. 14 and 4. 15 ttiat at this higher Reynolds nuoiheri 
tangential velocity and temperature profiles are almost 
identical for the present computations and the VSL model. 
It is due to the fact that at higher Reynolds numher the 
terms neglected in VSL model contribute very little. 
However, Figure 4. 16 shows that even at this Reynolds 
number, pressure profiles are different in the vicinity of 
the wall. This is due to the viscous terms dropped in case 
of the VSL model in the noinnal momentum equation. The 
difference in the density profile near the wall is not 
visible due to the scale (Figure 4.17). At this Reynolds 
number, post-shocK pressure, ten^erature and tangential 
velocity and shock stand-off distance are the same from 
the two approaches. 

Comparison with Experimental and other Theoretical H^esults 
Computations were carried out for the VSL model and 
the present method in order to compare them with available 
experimental and other theoretical results. Russel (1966) 
measured the density profiles in a rarefied stream. Figure 
4. 16 shows the comparison of density profiles from these 
two methods with the experimental data of Russel for the 
cold wall case ( Tw = 0. 26 To). It is noticed that the 
density profile computed from the full Havier-Stokes 
equations is in much better agreement with experimental 
data of Russel tlian tHe profile from VSL equations. Figure 
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4.19 shows density profile for the case of hot wall ( 
Tw s To). Once again it is seen that the density profile 
as predicted by the Navier-Stokes equations is much closer 
to the experimental results. 

Figures 4. SO and 4. Si show the comparison of density 
and temperature profiles with Honte Carlo results of 
Vogenltz and Takata (1970). For this analysis, Prandtl 
number of the gas has been taken to be S/3 which is 
consistent with the hard sphere model used for the Honte 
Carlo simulation. Figure 4. SO shows the comparison of the 
density profiles. The figure also Includes the results of 
Levlnsky and Yoshlhara (196S). It Is seen that the present 
computations agree well with Monte Carlo results while 
Levlnsky and Yoshlhara underpredict the density in 
comparison to the Honte Carlo results. However, Figure 
4. SI shows that resent results are closer to results of 
Levlnsky and Yoshlhara and both lu'edict higher temperature 
in comparison to Honte Carlo results. 

Eao (1964b) computed the hypersonic flow in the 
stagnation region of a blunt body using the concept of 
local similarity. Figures 4. S8 and 4. S3 show the 
comparison of profiles for the case studied by Cao. It is 
seen from Figure 4. SS that the present computations 
predict smaller value of shock stand-off distance in 
comparison to results df Hao. This figure also shows the 
hypersonic merged layer profile (^ain and Prabha 1963). It 
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IS seen tHat though, the normal velocity profiles are 
reasonably close to one another from th*'ee approaches, the 
behaviour of tangential velocity is quite different in 
case of the merged layer solution. It is due to the fact 
that in case of merged layer, the boundary condition is 
being applied at the frees tream while in other two methods 
it is being applied at the shock. In case of tangential 
velocity, the boundary condition is same in all the three 
methods, namely u = i, while for normal velocity, 
temperature et cetera, it is different in these methods. 
Figure 4. 23 shows the variation of temperature and 
pressure in the flow field. It is seen that temperature 
profile agrees well from these three methods, but pressure 
profiles differ. This difference in the pressure profile 
la because Eao has used inviscld pressure distribution to 
evaluate second order derivative of pressure for 
integrating the continuity equation. In case of merged 
layer solution, this difference in the pressure profile is 
due to the effect of shock transition zone. It seems that 
the Influence of the shock transition zone is more 
pronounced on pressure than on temperature or normal 
velocity. 

Computations were performed for the conditions used 
by Boy 1 an <197l) in the experinmnial measurement of heat 
transfer coefficient. Figure 4. 24 showg the variation of 
heat transfer coefficient with parameter e s i/ifRe and 
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conpax^lson witli cxperlncnlai data of Boylan. VSL pasult-s 
for tlie same conditions are also presented. It is seen 
tliat present neat transfer predictions are closer to tlie 
experimental data even for lower Reynolds number. While 
VSL results deviate from the experimental data in the 
lower Reynolds number regime. This again shows that the 
full Havler-Stokes equation model is better than VSL model 
for heat transfer predictions. 

Figure 4. 25 shows the influence of rarefaction on the 
heat transfer coefficient. In this case, computations were 
performed for = 20. The data is compared against the 
experimental data of Whlttllff and Wilson, as given by 
Vidal and Whlttllff (1963). It is seen that for rarefaction 
parameter as low as 0. 1, the present computations yield 
results which are quite close to the experimental data. 
Figure 4.26 shows the similar comparison against more 
recent experimental data of Coleman et al. (1977). In this 
case also it is seen that the present computations give 
good results for rarefaction parameter even below o. l. 

Figure 4. 27 shows the variation of the recovery 
factor with post shock Reynolds number. This has been 
compared with the experinmntal data of Hickman and Qiedt 
(1963). This figure also shows the recovery factor as 
conyputed from VSL model. It is seen that the recovery 
factor computed from present model reflects the actual 
behaviour quite accurately while the recovery factor 
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computed from VSL model does not five even tne qualitative 
behaviour correctly. 

Influence ^ Surface Temperature on Heat uitii? 

Computations were performed to estimate the heat 
transfer coefficient using the present approach for = 

5, 10, 16 and 20 and for wall temperature ratio Tw/To s 
0. 2, 0. 4, 0. 6 and O. 8. Figure 4. 28 shows the variation of 
heat transfer coefficient as function of Cheng's 
rarefaction parameter figure shows that, for a 

given surface temperature, Hach number variation from 6 to 
20 has negligible effect on the heat transfer coefficient, 
rmis is an expected behaviour in view of the Hach number 
independence principle applicable for hypersonic flow. 
However, it is noticed that there are distinct curves for 
different wall temperatures. These curves tend to collapse 
into single curve at higher values of From this It 

can be concluded that Qieng's rarefaction parameter falls 
to completely correlate the heat transfer behaviour as the 
level of rarefaction increases. Temperature of the surface 
also needs to be taken into account as already noted by 
Jain and Eumar (1986). This is despite the fact that 
rarefaction parameter u a function of wall 

temperature through the viscosity and Chapman-Rublsen 
constant. 

Usually aurface slip results in redaction of heat 
flux to the surface. Oavis <1970) has shown the influence 
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Of botli surface slip and stiocK slip on tlie beat transfer 
coefficient using VSL equations. Both of these effects 
tend to decrease the beat flux to tbe surface. Jain 'and 
Xdlmurtby (1974), fdio Used a full Havler-StoKes noerged 
layer model reported tbe increase in beat flux due to 
surface slip, in order to find out wbetber tbis phenomenon 
of increasing beat flux is associated with Havier-stoKes 
equations, computations were performed for tbe free stream 
and surface conditions used by Davis. It was found that 
influence of surface slip was to decrease tbe beat flux 
even In tbe present model of Havler>StoKes solution. In 
order to resolve this problem, computations were performed 

for different surface temperatures because surface 

\ 

temperature conditions were different In case of Jain and 
Adlmurtby. Figure 4. 29 shows tbe variation of beat 
transfer coefficient with and without surface slip for two 
different Reynolds numbers against wall temperature. It Is 
seen that surface slip increases tbe beat flux for lower 
surface temperatures. But at relatively higher surface 
temperatures, surface slip leads to reduction in beat 
flux. Surface slip reduces tbe temperature gradient at tbe 
wall (Figure 4.30), but at tbe same tloie increases tbe 
temperature at tbe base of shock layer leading to 
increased thermal conductivity. At lower surface 
temperatures this increase in thermal conductivity in 
combination with tbe redaction in temperature gradient 



142 


results in an increase in tne heat flux. At higher surface 
tenperatures this increase in thermal conductivity is not 
sufficient to compensate for the reduction in temperature 
gradient. 

4 . 3 Kesults for Chemical iv Heactins Gas 

Computations were performed for a chemically reacting 
gas with variable specific heat. The free stream 
conditions and vehicle parameters are representative of 
Shuttle orbiter nose during re-entry. These conditions and 
parameters are given In Table 4. The confuted results are 
compared with the merged layer results ( Gupta and 
Slmmonds 1986) and DSHC results < Moss and Bird 1985). In 
these computations Prandtl number and Lewis number were 
taken to be constant as in case of the merged layer 
analysis. The confutations were performed for non- 
catalytlc wall case which Is more realistic physical 
situation for Shuttle. 

Figure 4. 31 shows the comparison of heat transfer 
coefficient variation with Cheng's rarefaction parameter. 
It Is seen that present results agree well with the merged 
layer results and DSHC results except for 109 km altitude 
case. This difference in the heat transfer rates between 
present fozmmiiatlon and merged layer formulation at higher 
altitude is possibly due to the neglect of tangential 
transport effects in the shock transition zone which is 




approximately toz of tue shock layer thickness even at 
92. 35 Km altitude according to the oterged layer 
computations. Barller, similar results < Jain and Kumar 
1967) were computed for these conditions for perfect gas. 
It was noticed that the computed results deviated from the 
DSHC results at an altitude of 104 Km. For present 
coimputations the agreement is good upto 109 Km. For 104 Km 

altitude, the flow is frozen and it was expected that 

\ 

perfect gas results will hold good. This deviation is 
caused hy variahle specific heat of the gas which is not 
accounted for in the perfect gas computations. 

Figure 4. 32 shows the variation of temperature across 
the shock layer for 92. 35 Km altitude case. The profile 
has heen compared with merged layer profile and DSHC 
profile. It is seen that within the shock layer, the 
profiles are in good agreement with one another. 

Figures 4. 33 through 4. 37 show the variation of mass 
concentration across the shock layer for various species. 
The profiles have heen compared with corresponding 
profiles for merged layer and DSHC results. It is seen 
that in general, present profiles are in better agreement 
with merged layer profile. 

Figure 4. 33 shows the variation of mass concentration 
of molecular nitrogen. It Is seen that at the edge of the 
shock layer the concentration la same for both the present 
method and the merged layer solution. There is good 
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agreement In the profile across the shocK layer between 
present predictions and merged layer predictions. The 
continnnm methods predict higher concentration of 
molecular nitrogen as compared to DSHC results. 

Figure 4. 34 shows the variation of mass concentration 
of molecular oxygen across the shock layer. Once again it 
is noticed that there is a good agreement in the profiles 
predicted hy present computations and merged layer 
solution wd these concentrations are higher than the 
concentrations predicted hy DSHC. One noticeable aspect of 
this flgxire is that at the wall the concentration of 
molecular oxygen is 0.025 for continuum methods while for 
DSHC it is virtually zero, mis implies t^at there is 
conplete dissociation of molecular oxygen at about 2000*^1:. 
This is because near the wall, the concentration is 
controlled by diffusion and not by chemical reaction. 

Figures 4. 35, 4. 35 and 4. 37 show the variation of 
mass concentration of atomic oxygen, atomic nitrogen and 
nitric oxide respectively. It is seen that there is good 
agreement between present results and merged layer results 
while DSHC predictions are higher than these 
concentrations. This difference in concentration levels of 
various species as predicted by continuum methods and DSKC 
method can be attributed to difference in dependence of 
Viscosity and chemical reaction rates on temperature as 
used in two different methods. It is also to be noted that 



DSHC results predict blflier concentrations of atojr,ic 
species and lower concentrations of moleculeu' species. 
Possibly It Is due to inadequate accounting of multi -body 
collision processes, wblcb are mainly responsible for 
recombination reactions. 

It is noticed from figures 4. 33 to 4. 37 tbat various 
concentrations from tbe two continuum methods, namely the 
shocK layer solution using full Havler-StoKes equations 
and merged layer solution are In good agreement . For this 
altitude, the merged layer computations show that shock 
transition zone is approximately 70'/. to 60Z of the shock 
layer thickness. Even though, this thickness is comparable 
to the shock layer thickness, still there is a good 
agreement in the post shock concentrations for the two 
methods. It shows that the effect of chemical reactions in 
the shock transition zone on the mass concentrations of 
various species is not very significant. 

Figures 4. 3d and 4. 39 show the temperature and mass 
concentration profiles for the free stream conditional 
corresponding to 104 Em altitude. It Is observed from 
Figure 4. 38, showing the temperature profile, that present 
method lo'edlcts higher temperature across the shock layer 
in conqparlson to the merged layer and bSHC results. This 
higher temperature prediction by present method is 
possibly due to shock conditions, which assume thin shock 
and thus neglect transfer of momentum and energy parallel 



to til© shock. This effect can he significant here as the 
shock transition zone is approximately 6 times the shock 
layer thickness according to the merged layer 

computations. Figure 4.39 shows the concentration profile 
of the molecular nitrogen and molecular oxygen. It is 
noticed that the concentrations are constant across the 
shock layer. This implies that flow is frozen at this 
altitude for the conditions used in the analysis. 
influence of Varlahle Transport Properties 

In the results presented so far, Prandtl number and 
Lewis number were assumed to be constant. In order to 
study the influence of the variation of Prandtl number and 
Lewis number on the behaviour of the flow parameters, 
computations were carried out for Shuttle conditions 
considering variable Prandtl niimber and Lewis number. For 
Lewis number estimation, binary diffusion of atomic oxygen 
into molecular nitrogen was considered as suggested by 
Koss (1971). Figure 4.40 shows the variation of heat flux 
with Cheng's parameter for different flight conditions. 
!nie mixture properties were evaluated using Armaly and 
Sutton (1980,1981) correlation. It is seen that the use of 
variable Lewis number In the shock layer as a function of 
temperature and composition does not affect the heat flux 
appreciably. On the other hand, use of variable Prandtl 
number in the shock layer has significant effect on the 
heat flux. For 92. 36 im altitude case the heat flux 
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increases by approximately 20*/ . Similar computations were 
carried out using WllKe's correlation also for tbe mixture 
properties. Figure 4. 41 shows the variation of heat flux 
with Cheng's parameter for this case. In this case it is 
noticed that the use of variable Prandtl number and Lewis 
number in the shocK layer has no appreciable effect on the 
heat £lux to the wall. Table 5 shows the heat flux values 
for the various cases considered in this study. In case 
Armaly and Sutton correlation is used, the increase in 
heat flux is associated with the increased value of 
thermal conductivity. 

Figures 4. 42 through 4. 49 show the variation of 
different, flow variables across the shocK layer for the 
two mixture property correlations. The profiles agree well 
for the two cases except for marginal quantitative changes 
in the values of the flow variables. Figure 4.42 shows the 
variation of the temperature across the shocK layer. It is 
seen that use of Armaly and Sutton correlation predicts 
higher teim;>erature at the wall and lower temperature at 
the ShocK in comparison to the corresponding predictions 
for WilKe's correlations. This is because Armaly and 
Sutton correlation predict higher values of thermal 
conductivity at higher temperature and thus causes more 
thermal diffusion in the post shocK region leading to 
lower temperature there. This heat transfer leads to 
higher temperature near the wall. 
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Plgxire 4.43 shows the concentration profile for 
molecular nitrogen. This figure shows that use of Aroaly 
and Sutton correlation predicts the higher concentration 
of molecular nitrogen in the shock layer the 
corresponding predictions for Wilke's correlation. This is 
due to lower temperature predictions in case of Armaly and 
Sutton correlation and thus lower rate of dissociation. 
Hear the wall, where the gas temperature is not high 
enough for dissociation to take place, and flow is 
controlled hy diffusion resulting in significant 
concentration of H and lower concentration of Hg. Figure 
4.44 shows the concentration profile for molecular oxygen. 
Here it is noticed that use of Armaly and Sutton and 
Wilke's correlation give similar values of molecular 
oxygen concentration. This is due to fact that in the 
entire shock layer except near the wall, the temperature 
is high enough for oxygen molecules to dissociate. Thus 
both the methods give similar values of concentration of O 
and 0£. Hear the wall, the concentration is governed by 
the diffusion. 

Figures 4. 45, 4. 46, 4. 47 and 4. 46 show the 
concentration profiles for atomic oxygen, atomic nitrogen, 
nitric oxide and ionised nitric oxide, it is seen that use 
of Wilke's correlation in the computations leads to 
prediction of higher concentration of H and HO'*’ than the 
corresponding predictions for Azmialy and Sutton 
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correlation, nils is tue consequence of lilgber temperature 
in the shocK layer for the case when WilKe's correlation 
IS used. 

Figure 4.49 shows the variation of Prandtl nuDOer in 
the shock layer. It is seen that use of Wilke's 
correlation .leads to higher values of Prandtl numher 
throughout the shock layer. This is mainly due to 
underprediction of thermal conductivity when Wilke's 
correlation is used and partly due to chazige in local 
concentration of various species. This change In thermal 
conductivity causes increase in heat flux to the surface 
when Armaly and Sutton correlation is used. Table 6 gives 
the comparison of wall values and post shock values of 
different flow variables for the two cases. 

Influence of Bach Bunber on Heat Transfer to Wall 

In hypersonic flow, according to Hach number 
Independence principle, the heat transfer coefficient is 
Independent of Mach number. It is true for the case of 
perfect gas as already observed in figure 4. 26. However, 
in case of chemically reacting gas flows, Hach number 
Independence principle may not hold. Higher free stream 
velocity implies higher post shock temperature and thus 
higher temperature in the shock layer. The reaction rates 
being dependent on the actual temperature, chemical 
reactions (dissociation, Ionisation et cetera) will 
proceed at different rates for different free stream 
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velocities and tbus result in different neat fliuc levels. 
Tnis variation of velocity will not result in any change 
In Reynolds number based on post snocK temperature or 
based on reference conditions used in the present worK if 
the square-root viscosity law is used. However, in actual 
case, Reynolds number based on reference conditions or 
based on post shock conditions will not be constant but 
change slightly with velocity as square-root viscosity law 
is not strictly valid. 

In order to study the Influence of free stream 
velocity on heat transfer coefficient, computations were 
carried out for Shuttle conditions with free stream 
velocity varying between 6 Cm/sec to iO xm/sec. It should 
be noted here that for higher free-stream velocity, the 
seven species air model may not be appropriate and better 
model may be necessary. Figure 4. 50 shows the variation of 
heat transfer coefficient against free stream velocity for 
different altitudes. It is seen that for altitude of 
104. 93 Xm and 109. T5 Km. where flow is frozen, heat 
transfer coefficient decreases with increase in free 
stream velocity. Heat transfer coefficient decreases with 
increase in free stream velocity for 92. 36 Km altitude 
while for altitude 99.49 Dn heat transfer coefficient 
increases with increase in free stream velocity. In case 
of 104. 93 Km and 109 Xm altitude, the flow is frozen and 
one might expect that heat transfer coefficient will 



remain cons-tant with free stream velocity as Hach number 
independence principle should apply in these cases. The 
variation in heat transfer coefficient is due to the 
temperature dependence of specific heat, ihe variable 
specific lie at leads to lower non-dimensional temperature 
in the sliock layer and thus lower heat transfer 
coefficient. Figure 4 . 51 shows the variation of heat 
transfer coefficient against Reynolds number based on 
reference conditions, it is seen that heat transfer 
coefficient varies very sharply with Reynolds number as 
there in not much change in Reynolds number with change in 
free stream velocity. This leads one to conclude that free 
stream velocity and thus the actual temperature in the 
shock layer should be accounted for while considering the 
Reynolds number parameter in case of chemically reacting 
flows. 

In order to estimate the range of the Reynolds number 
for which present scheme converges for chemically reacting 
flow, computations were carried out for AFE trajectory. 
The trajectory details are given in Table 7. It should be 
noted here that at these high free-stream velocities, 
seven species model for the air is not likely to 
represent the actual situation for non-equilibrium flow. 
Moreover radiation is likely to influence the flow-feild 
significantly. severtheless, the computations were 
performed to assess the capability of the present scheme. 
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Figure 4.52 shows the variation of stagnation point heft 
transfer coefficient with post shocK Reynolds niimber for 
APE vehicle. The computations converged without any 
problem for the entire trajectory range, which corresponds 
to 

3 < Res < 1760 

This upper limit of Re^ 1760 is achieved for the present 
scheme as against of Re, 600 for VSL scheme (Moss 1971). 
The details regarding the convergence of various surface 
quantities and profiles for higher Reynolds number case 
have already been discussed in Chapter 3. 

Figure 4. 52 shows the variation of stagnation point 
heat transfer coefficient with post shocK Reynolds number 
for AFB vehicle. This figure shows the heat transfer 
coefficient for frozen flow and both for catalytic and 
non-catalytic wall in case of non-equili)>x*ium flow. It is 
observed that heat flux is lowest in case of non-catalytic 
wall idxile it is highest for catalytic wall. In case of 
non-equilibrium flow, higher temperature induces 
dissociation and ionisation of various species. This 
process absorbs heat and lowers the temperature of the gas 
in the shock layer. This lower tea4>erature results in 
lower heat flux, however, in case of catalytic wall, ell 
the ionised and dissociated species combine to give 
molecular species near the wall and release the energy 
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associated nflth. dissociation and ionisation. This leads to 
increased heat flux to the wall. 

Figure 4. 53 shows the variation of wall pressure with 
post shock Reynolds number for APE vehicle. It is observed 
that non-equilibrium flow predicts higher wall pressure in 
comparison to frozen flow. Due to chemical reactions, 
various species dissociate and ionise. This reduces the 
molecular weight of the gas resulting in increased 
pressure in case of non-equi 1 ibrium flow. 

It was stated earlier that present computations 
converged for AFE vehicle trajectory at an altitude of 
78, 9 Km corresponding to post shock Reynolds number of 
1780. As this Reynolds number is reasonably high, it was 
expected that flow may be in chemical equilibrium. In 
order to check whether or not present computations capture 
the equilibrium composition, the same was estimated 
corresponding to the pressure and temperattire profiles 
obtained from the non-equi librium computations. Figures 
4. 54 to 4. 59 show the concentration profiles of Hg, Og, 
HO, O, H and HO* respectively as obtained by non- 
equilibrium computations and their comparison with 
corresponding equilibrium concentration for the pressure 
and temperature of the non-equilibrium flow. It is 
observed that in most of the shock layer region, 

equilibrium and non- equilibrium concentrations are 
identical. These concentrations differ near the wall and 



ncax* me sliocK* At tlie wall, non- e<iui librium computations 
give liiglier values of atomic species and lower values of 
molecular species concentrations in comparison to 
equilibrium values. Tbis is due to the fact that in case 
of non-equi 1 Ibrium computations, concentrations near the 
wall are controlled by diffusion while in case of 
equilibrium, concentrations are controlled by chemical 
Kinetics. Hear the shocK, the difference in concentrations 
is due to the boundary conditions to be satisfied by the 
concentrations at the shocK. The difference in the 
concentration profile near the wall contributes to the 
heat flux to wall. It is due to energy transport 
associated with diffusion of species. It implies that even 
if the flow is mostly in chemical equilibrium and non- 
equilibrium computations do capture equilibrium 
concentration profiles as 'seen in figures 4. 54 to 4. 59, 
the equilibrium and non- equilibrium heat flux to the wall 
will differ. This difference in the heat flux will grow 
with increase in the Reynolds number as it will result in 
sharper teo^erature gradient and therefore concentration 
gradient near the wall. Moreover, when equilibrium code is 
run, this diffrence in concentrations near the wall and 
the shock will influence flow field in the entire shock 
layer and the heat flux to the wall. 
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Free Stream Conditions and Yeiilcle Parameter 


Free Stream conditions 


Alt. 

Cm 

DensltF 

Cg/m3 

u®* 

Cm/sec 

T®* 

Or 

Hole 

XH2 

Fractions 

X02 

XO 

92. 35 

2. 164e-6 

7. 50 

180. 0 

. 783 

. 217 

000 

99.49 

5. 906e-7 

7. 50 

190. 0 

. 783 

. 217 

000 

104. 93 

2. 457e-7 

7.47 

223. 0 

. 782 

. 153 

065 

109. 75 

1. 146e-7 

7. 47 

249.0 

. 771 

. 123 

106 

Yeblcle Farameter wltli Altitude 

Altitude 

Hose Radius Rl> 

wall 

Temperature 

cm 


m 



Oc 


92. 35 


1. 296 


1043. 0 


99.49 


1. 362 



800. 0 


104. 93 


1. 362 



560. 0 



109. 75 


1. 362 


4E0. 0 





Indllnence of Property Variation on Heat Flnr 


(a) Uae of Amalp and Sutton Correlation 
Heat Flux in rw/m2 


Pr 

Le 

Altitude (Km) 

92. 35 99. 49 104. 39 

109. 0 

0. 7 

1. 4 

92. 92 

47. 85 

30. 08 

19. 20 

0. 7 

Le(T) 

92. 10 

47. 64 

30. 08 

19. 21 

Pr(T) 

1.4 

110. 00 

51. 00 

32. 80 

19. 70 

Pr(T) 

Le(T) 

109. 26 

51.46 

32. 80 

19. 73 


(b) Hae of VllKe*s Correlation 




Heat Flux in Kw/id2 





Altitude (im) 


Pr 

Le 

92. 35 

99. 49 

104. 39 

109. 0 

0. 7 

1.4 

90. 20 

47. 10 

29. 86 

19. 10 

0. 7 

Le(T) 

89. 80 

47. 17 

29. 86 

19. 10 

Pr(T) 

1. 4 

90. 66 

47. 20 

30. 16 

19. 11 

Pr(T) 

Le(T) 

90. 20 

47.20 

30. 20 

19. 12 
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’TASIaEB 6 


Wall and PoBt snock Valnea o| Yar;al>lea 
AwaalY J Sutton and Wilke* a Correlation 



Correlation 

Armaly a Sutton 

Wilke 

Wall Valnea 

Temperature (^) 

2386. 0 

2043. 0 

Che 

0. 6604 

0. 6780 

Co2 

0. 0851 

0. 0243 

Cho 

0. 0553 

0. 0624 

Cq 

0. 1859 

0. 1880 

Ch 

0. 0532 

0. 0569 

Cho+ 

0. 000126 

0. 00013 

Pr 

0. 6711 

0. 6876 

Post Shock Valuea 

Temperature (9c) 

15999.0 

16267. 0 

Hormal Velocity 

im/sec) 762. 67 

768. 35 


0. 7153 

0. 7133 

Co2 

0. 1536 

0. 1539 

Cho 

0. 0216 

0. 0213 

Co 

0. 0754 

0. 0752 

Ch 

0. 0341 

0. 0361 

ChO+ 

0. 000068 

0. 000064 

Pr 

0. 6210 

0. 7200 
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TRAJBCT08T DETAILS FOR APB VBBICLB 


Alt. 

Bm 

Density 

X8/IS3 

UflO* 

Zm/sec 

Ok 

Hole Fractions 

XH2 X02 XO 

il4. 3 

6. 306e-8 

9. 905 

285. 8 

. 7571 

. toil 

. 1416 

liO. 5 

9. 669e-8 

9. 909 

251. 2 

. 7695 

. 1205 

. 1100 

107. 2 

1. 661e-7 

9. 911 

230. 8 

. 7773 

. 1405 

. 0622 

103. 9 

2. 876e-7 

9. 913 

214.4 

. 7832 

. 1581 

.0587 

100. 6 

5. 309e-7 

9. 915 

203. 7 

. 7848 

. 1740 

.0412 

97. 7 

9. 056e-7 

9. 912 

196. 8 

. 7865 

. 1864 

. 0271 

94. 8 

1. 513e-6 

9. 909 

191. 7 

. 7870 

. 1977 

.0153 

92. 3 

2. 392e-6 

9. 898 

188.4 

. 7873 

. 2049 

.0078 

89. 9 

3. 712e-6 

9.886 

186. 7 

. 7880 

. 2120 

. 0000 

87. 7 

5. 350e-6 

9. 863 

189. 5 

. 7880 

. 2120 

. 0000 

85. 6 

7. 615e-6 

9. 832 

191. 5 

. 7880 

. 2120 

.0000 

83. 9 

1. 015e-5 

9. 785 

193. 0 

. 7880 

. 2120 

.0000 

82. 4 

1. 364e-5 

9. 539 

194.0 

. 7880 

. 2120 

. 0000 

81. 2 

1. 614e-5 

9. 653 

194.4 

. 7880 

. 2120 

. 0000 

80. 1 

1. 927e-5 

9. 569 

196. 5 

. 7880 

. 2120 

. 0000 

79.4 

2. i88e-5 

9. 471 

195. 6 

. 7880 

. 2120 

.0000 

78. 9 

2. 388e>5 

9. 365 

196. 2 

.7880 

. 2120 

.0000 


Wall Temperature { till 

Bose Radius : 2. 53 m 
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Fm 4.2 Shock Stand-oft Distance reft Hot Wall Case 
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Data or RussTl for Cold Wall 
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CHAPTEiR S 
CONCLUSIONS 


5* * IntrcductjLop 

In this ctiapter, the results of the present 
investigation are summarised. These results have already 
been discussed in Chapter 4, where the detailed results of 
tbe present worK were presented and were con«>ared with 
other continuum, DSHC and experimental results. The 
salient features of these results are reiterated here with 
a view to evaluate how far the questions raised in section 
i.4 have been answered. These questions are; 

(i) IS It possible to estend tbe ranse o* Tslldlty of 
sbocK layer metbod by use of full Havler-stoKes equations? 
<ii) HOW do tbe cbemical reactions witbin tbe sbocK 
transition zone influence tbe sbocR layer flow-field at 


low Reynolds nuslbers, where tbe sbocK ware is tblcR? 

,lli) IS it possible to extend tbe r»ae of Reynolds number 


for wblcb non-equlllbriuro computations converfe? 

(iv, HOW to overcome tbe violation of mass conservation 
law nomally occurin, in numerical studies of ob«.lc.lly 


reacting flows? 
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(V) Is it possible to oyercoi&e tbe nunerioel instalsitity 

✓ 

caused by spurious negative concentrations? 

The present work is ootivated by the fact that most 
of tHe continuum analyses for the chemically reacting flow 
are being carried out by using either VSL/PHS method -or 
merged layer approach. The governing Havi«r-*stoKes 

equations are simplified in the VSL or PHS method- Certain 
terms are dropped from Havier-StoKes. equations to '^simplify 
the governing equations. These terns are not negligible-! at 
low Reynolds number range and thus these methods fail to 
give the correct behaviour of flow structure and surface 
characteristics at low Reynolds nuntber. This llmlirs the 
applicability of the VSL and PHS methods. The merged lairer 
method requires the solution of the Havler>-St<4ces 

equations from the body surface to the free-stream through 
the shock, which comes out as the part of the sofutidn. 
The difficulties associated with the numerical integration 
of the Havier-Stokes equations make the numerical 
computations very expensive. Horeover, at higher Rgrynoids 
numbers, the shock becomes very thin making it difVlrdU.lt 
to resolve the shock structure unless a very fine gi^id ^^Is 
chosen. This limits the applicability of the merged ' layer 
^thod. Moreover, the numerical problenm associated - wdth 
the dhemiodily; reacting flow, namely the nuadridal 
1 violation of overall mass conservation due to 'htiidiKlhal 
, reactions and instability caused by spubious i negative 
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concentrations, also need to be resolved. 

In tne present work, fun Havier-Stokes equations for 
cliemicallY reacting flow are integrated in tne stagnation 
zone. THe equations are integrated numerically between tbe 
body and tbe shock using sm^face slip and the shock slip 
boundary conditions. ^Hie in^ortant results are briefly 
discussed in the following sections. 

5. 2 Extension Validity of Shock Laver Approach with 
Full Mavier-Stokes Bquations 

It is known that the stagnation point heat transfer 
rates predicted by Cheng's two thin layer theory agree 
well with the experimental values even at very low 
Reynolds nuiifl>ers. However, experimental investigations by 
Ahouse and Bogdonoff (1969) on detailed flow field in the 
stagnation region. Honte Carlo simulation of the flow 
field at low Reynolds nuisber( Vogenitz and Takata, 1970} 
raised the doubt on the validity of the two thin layer 
approach. Ahouse and Bogdonoff (1969) attributed this 
discrepancy to the failure of continuum approach. This 
failure was attributed by. Jain and Adimorthy (1974) to the 
thin layer approximations. The present solution, using 
full ; |iavi4^'>iStokes equations in the shock layer, lOiovs a 
fairly good agreement with the experimental data and Honte 
Carlo results at discussed in Chapter 4. The use of full 
Kavier-Stokes equations in the shock layer predicts 
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following additional features wUlcH the VSL model fails to 
predict even qualitatively . 

(i) Behaviour of the predicted wall inr^esstire agrees well 
with the experimentally observed behaviour of wall 
pressure at low Reynolds number. 

(ii) Theoretically predicted recovery factor agrees well 
with the experimentally observed recovery factor in the 
low Reynolds number regime. 

(iii) There is a reasonably good agreement between the 
predicted temperature and density profiles and 
experimental data. 

(iv) It is also observed that at higher Reynolds number, 
the computed heat transfer coefficient agrees well with 
boundary layer prediction for the heat transfer 
coefficient. However, this phenomenon is also true of 
shock layer model. 

These factors show that the use of full Havier-StoKes 
equations reproduces many ln®»ortant features of the 
hypersonic rarefied flow. At the same time at higher 
Reynolds nuna>er, this formulation yields the results which 
are in good agreement with boundary layer theory 
predictions. 

'f>'? anfmeiiye CheleaJ pegctlgai lai ^fi£l 

W9m m mt iMM 

in t|#<Mipreient analysis, derivation of shock slip 

; 1 * ' ' ' 

condltlons/post shock conditions assuans an infinitely 
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sliocK. nxis tacitly assiunas tliat clxemlcal reactions 
in the shock transition zone do not influence the post 
shock chemical composition. The thickness of the shock 
transition zone Increases as the Reynolds number 
decreases. Due to high temperature in the shock transition 
zone, the reactions are f^st. With increased thickness of 
the shock transition zone, it is possible that post shock 
composition will be different from the free stream 
composition. 

Merged layer and DSHC methods solve the flow field 
through the shock transition zone and thus Include the 
Influence of chemical reactions in the shock transition 
zone on the flow field. It is observed that the computed 
profiles agree well with the profiles obtained form merged 
layer analysis and DSKC results. This shows that the 
Influence of chemical reactions in the shock transition 
zone on the flow field is not significant. The influence 
of chemical reactions in the shock transition zone, if 

any, is limited to the neighbourhood of the shock. 

* , 

S. 4 Copvereence »on-Bqqi I ibr ium Computational Procedure 

The computer code based on present fomulation was 
rdh for a wide range of free-stream conditions. So 
convergence prohlen was noticed for AfV iraiectory for 

3 i tag s iT6d 


t 


2£3 

wliere is tlie Reynolds nuobez* ]>ased on tlxe post. slxocK 
conditions. It is observed tbat non'*e<iui 1 ibrium 
computations capture equilibrium concentrations in most of 
tbe sbocK layer at bigber Reynolds number. The upper limit 
of Rcj for convergence of non-equilibrium code for viscous 
SbocK layer (VSLy model was reported to be 500 by 
lloss(i97i}. Thus it is seen tbat use of full Havier-StoKes 
equations along witb UDR scbeix» and present treatment of 
mass production terms bas increased tbe range of Reynolds 
number for wbicb non- equilibrium analysis converges. 

5. 5 Hass Conservation and Species Concentrations 

Tbe violation of mass conservation noraally oceuring 
in tbe numerical studies of cbemically reacting flows, is 
taKen care of by considering the elemental mass 
conservation at each iteration as explained in section 
3. 4. This treatment of production terms ensures tbat at 
each iteration, overall mass is conserved vitb change in 
concentration of various species. UlIs treatment also 
ensures tbat overall laass is conserved with change in 
temperature in tbe presence of chemical reactions. 

At abf iteration, occurrence of spurious negative 
cbboen^rations will result in improper estimation of 
;par*^|i;fi|eti4i|. ..telm This may lead to breakdown of solution 
prooednre, ^sbiis is taken care of by speciffUni tbe beun^ 
for tbe concentration of <tliie vnrions ipooioi^ 'Hm 
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concentrations are forced to be positive at any iteration. 
An upper limit lias been specified for ail tliese 
concentrations. Tbis upper limit is estimated from tbe 
free stream composition of tbe air. 

S. 6 Some froblems for Furtber RcaearcH 

The success of the full Havier-Stokes formulation in 
the shock layer method in describing the chemically 
reacting hypersonic flow in the stagnation region of a 
blunt body, encourages the following extensions to the 
present work. 

(i) Inclusion of radiative heat transfer effects in the 
present foxmnilatlon for blunt body flow problem. 

(ii) Introduction of mass injection 

(iiil Introduction of multi •sompoaent dif^asica and the 
influence of thermal and pressure diffusion. 

(iv) Introduction of chemical reaction between gas and 
surface material. / 

(V) Investigation of influence of higher order 
approximations on the stagnation region. 

(vl)!Che present solution is limited to the stagnation 
region of a blunt body. This analysis can be extended to 
4 nt*ia*ate aawler-Stokes e<iuations in the downstream region 

of the blunt boidy with present solution taken as initial 

, '( 

Si . . , 
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AFPBNDIX 1 


IBFLUEHCB LIBKAPT^ATIOI Of PKOPOCTIOH TBHH 


Ttxougli the analysis involves multi -component multi- 
reaction gas mixture, for the saKe of developing an 
understanding of the influence of linearisation of 
production terms, consider the case of single reaction 
given by 


Og 


Kb 


> 20 


( A. i. 1) 


For this reaction u<ass ^oduction rata is given by 


r C C -I 

£ H I p K -p® 3C ( — ^ I 

O L ^ Hop ^ Ho 


Hoa " Ho 

( A. 1. £ ) 


0 £ 


- C C n 

= - H r P ^ 

. o£ L ^ Hoa ^ 

< A. 1. 3 1 


Using the fact that Hoa = 2 Hq * equations ( A. i» i y 
Sfkd ( A. 1.2 y can be written as: 

wq = a t Coa “ o ^ 


( A. 1 . 4 y 
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and 


iro2 = A t a CqS - C02 J 


wnere 


( A. i. 5 ) 


A = 9 Xf 

« = 9 Xi^ / Xf 

At tills Stage, the s\ub of the two mass production 
rates is zero, that is: 


Wo + Woa s 0 


< A. 1. 6 ) 

whereas, linearisation of the equations ( A. 1,4 
( A. i. 5 ) yields 


and 


Wq = A Cpg - A a Cp 

( A. i. 7 ) 

Wpg = A a CqS - A Cpg 

( A, t.d ) 

where the underlined terms have taken their walues at the 
previous itextAtlon* 

Addition of the two production rate terms gives 
Wp ♦ Wpg s A < Cpg - Cpg )+A ) 

( A. 1.9 » 

Qm applfl>>Me mass conservation condition 
‘ i cp ♦ Cpg a I 


( A. I. 10 ) 
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one gets 

Co - ^ = - < Co£ - Cog) 

( A. 1. 11 ) 

Finally, equations ( A. l. ii ) and ( a. 1.9 ) yield 

Wo + ^02 = A < 1 + a ^ ) < Cpg - Cog > 

( A. 1. 12 ) 

WHICH is non-zero quantity until new and old values of 
Co 2 are Identical. 

a J Hows tHat 1 Inearlsatlon of production terms eatiai»« 
net S^ss production rate to jbe non- zero, this fact sHould 
toe taken into account wHlle t^ details g£ ttoe solution 
procedure are toeing worked out. 


